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PREFACE 


This  Study  is  a  continuation  of  the  work  done  by  previous 
researchers,  at  A.F.I.T,  to  investigate  the  heat  transfer  to  a  flat 
plate.  Much  of  the  focus,  however,  was  deliberated  to  the 
calibration  of  the  heat  flux  gages  and  the  effect  of  film  cooling 
on  the  rate  of  heat  transfer  to  a  flat  plate  in  a  shock  tube. 

Keeping  in  view  the  vicissitudes  involved  in  an 
experimental  project  of  this  nature  it  was  necessary  to  seek  help 
from  various  sources  and  people.  A  timely  abet  from  all  those  asked 
for,  indeed,  helped  me  in  finishing  this  project. 

My  greatest  gratitude  certainly  would  go  to  my  advisor  Dr. 
William  C.  Elrod,  who  at  all  times  was  available  to  help  me  at  some 
very  crucial  moments.  His  careful  supervision  not  for  once  allowed 
me  to  go  astray. 

It  had  mostly  been  the  efforts  of  Lt .  Col  Paul  I.  king  and  Mr 
Jay  Anderson  whose  experimental  expertise  got  me  through  the 
cumbersome  job  of  calibration.  I  owe  acknowledgment  to  them  for 
their  expertise  and  patient  listening. 

To  the  Laboratory  staff:  Mr  Nicholas  Yardich,  Mr  Andy  Pitts, 
Mr  Dan  Rioux,  Mr  Mark  Derriso;  I  am  extremely  grateful  for  the 
consistent  help  and  consultation  on  the  instrumentation  and 
hardware  throughout  my  research  study. 

Last,  but  not  the  least,  I  extend  my  thankfulness  to  my  wife 
Saeeda  and  my  lovely  daughter  Qurat-ul-Ain  for  their  equanimity. 
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ABSTRACT 


Owing  to  the  complexities  involved  in  studying  heat  transfer 
to  complex  geometries,  initial  research  is  carried  out  on  simple 
models.  The  results  of  such  initial  studies  can  be  used  to 
understand  some  complex  flow  patterns. 

During  this  study  a  simple  model  of  a  flat  plate  was  chosen  to 
determine  heat  flux  from  high  speed  flow  behind  a  shock  wave  into 
the  plate.  An  electrical  analog  circuit  was  used  to  convert  the 
output  of  the  heat  flux  gages  into  heat  flux. 

The  free  stream  turbulence  withen  the  flow  behind  a  shock  wave 
and  that  generated  by  a  turbulence  generator  was  determined  and  its 
effects  on  the  heat  transfer  data  were  then  studied  for  comparison 
with  previous  experimental  work.  The  main  focus,  however,  was  laid 
upon  the  calibration  of  heat  flux  gages  and  the  effects  of  film 
cooling  on  heat  flux. 

Thin  film  heat  flux  gages  were  calibrated  for  their  material 
properties.  The  calibration  constants  were  then  used  to  determined 
the  true  heat  flux  into  the  flat  plate. 

Film  cooling  flows  were  studied  under  three  different  primary 
flow  conditions.  Cool  air  was  injected  into  the  boundary  layer  on 
the  flat  plate  and  analyzed  for  its  effects  on  heat  transfer. 

Owing  to  the  transient  nature  of  the  flow  behind  the  shock 
wave  some  critical  properties  like  the  adiabatic  wall  temperature 
of  the  flat  plate,  and  the  heat  transfer  coefficient  could  not  be 
determined  during  the  film-cooling  studies.  As  a  result  this  part 


of  the  research  was  restricted  to  the  analysis  of  the 
experimentally  obtained  data. 


XVI 


TRANSIENT  HEAT  TRANSFER  MEASUREMENTS 


ON  A  FILM  COOLED  FLAT  PLATE 
IN  A  SHOCK  TUBE 

I .  INTRODUCTION 

Current  aviation  applications  demand  extensive  research  in  the 
field  of  gas  dynamics.  With  the  increasing  interest  of  gas  turbine 
designers  to  improve  upon  the  turbine  blade  technology;  researchers 
are  getting  more  and  more  involved  in  exploring  and  understanding 
the  basic  heat  transfer  processes  taking  place  in  this  extremely 
complex  zone  of  highly  turbulent  flow  region.  However,  in  order  to 
minimize  the  complexities  involved  in  the  heat  transfer 
measurements  of  turbine  cascade  flow  geometry,  researchers  have 
simplified  their  initial  stages  of  research  by  confining  themselves 
to  a  comparatively  facile  problem  of  the  flat  plate.  The  results  of 
such  a  study  can  lead  to  a  better  understanding  of  flow  over 
complex  geometries  like  the  turbine  blade. 

An  impulsive  rise  in  gas  temperature  can  be  achieved  by 
generating  a  normal  shock  in  a  transient  facility  like  the  shock 
tube.  The  measurement  of  heat  transfer  is  easier  in  a  shock  tube 
than  a  continuous  facility  like  a  hot  tunnel.  Models  in  hot  tunnels 
require  internal  cooling  to  establish  a  temperature  gradient  from 
which  heat  transfer  can  be  deduced.  It  must  also  be  capable  of 
withstanding  recovery  temperatures  continuously,  and  if  the  model 
has  a  large  thermal  capacity  the  tunnel  itself  must  be  operated  for 
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an  extended  period  of  time  so  as  to  achieve  uniform  conditions 
within  the  model.  In  contrast,  a  model  in  a  transient  facility, 
like  the  shock  tube,  never  attains  an  equilibrium  that  enables  the 
heat  flux  to  be  determined  (Shultz  and  Jones,  1973:3). 

Owing  to  the  short  duration  of  tests  in  the  shock  tube,  fast 
response  measurement  devices  are  required  to  collect  the  requisite 
information.  Such  a  collection  of  data  can  be  achieved  by  using 
fast  response  thin  film  heat  flux  gages  which  allow  transient 
affects  to  be  properly  sampled  and  recorded. 

An  increase  in  the  rate  of  heat  transfer  can  be  achieved  by 
generating  turbulence  within  the  flow.  However,  to  protect  a 
surface  exposed  to  high  temperature  environments,  an  effective  way 
of  reducing  heat  flux  to  the  operating  surface  is  required,  one 
such  way  is  by  using  film-cooling  techniques  which  not  only  reduce 
the  heat  flux  but  also  increase  the  component  life. 

BACKGROUND 

The  history  of  the  shock  tube  goes  back  to  the  1900' s. 
However,  its  widespread  use  has  been  facilitated  by  the 
technological  boom  during  World  War  ll.  The  slow  growth  of  the 
shock  tube  use  during  its  early  years  is  mostly  attributed  to  the 
lack  of  fast  response  instrumentation.  With  the  advancements  of 
piezoelectric  gages,  high  speed  photography,  and  an  increased 
research  into  supersonic  flow  phenomena  in  the  post  World  War  II 
era,  the  shock  tube  has  become  an  important  research  tool, 
especially  in  the  fields  of  chemical  kinetics,  dissociation  and 
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ionization,  aerodynamics,  and  heat  transfer  (Bradly,  1962; 
Hartunian,  et  al,  1959;  Gaydon,  1963). 

Lam  and  Crocco  (1959)  studied  the  characteristics  of  the  shock 
induced  unsteady  laminar  boundary  layer  on  a  semi  infinite  flat 
plate.  They  experimentally  found  the  existence  of  two  distinct 
domains:  one  near  the  shock  where  the  flow  is  quasi -steady  and  the 
other  an  unsteady  region  where  the  flow  characteristics  approach  a 
steady  state  asymptotically. 

Mirels  (1961)  presented  numerical  solutions  for  the  laminar 
boundary  layer  behind  a  strong  shock  advancing  into  stationary  air. 
Boundary  layer  transition  and  the  steady  rate  of  heat  transfer  to 
a  shock  induced  turbulent  boundary  layer  was  studied  by  Spence 
(1960)  and  Hartunian,  Russo  and  Marrone  (i960).  They  obtained 
correlations,  based  on  empirical  data,  for  the  rate  of  heat 
transfer . 

Davies  and  Bernstein  (1969)  investigated  the  shock  induced 
boundary  layer  on  a  semi  -  inf  ini te  flat  plate.  They  found  that  for 
the  laminar  boundary  layer,  with  zero  pressure  gradient,  the  flow 
is  steady  at  a  distance  x  from  the  leading  edge  of  the  flat  plate 
after  some  time  t  and  is  given  by  : 

X 

0.317 

where  U  is  the  velocity  of  the  flow  ,  outside  the  boundary  layer, 
behind  the  shock  wave.  For  a  turbulent  boundary  layer  they 
discovered  that  half  the  boundary  layer  length  is  steady. 

With  the  advancements  in  the  development  of  fast  response 
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instrumentation,  the  use  of  shock  tubes  has  further  been  enhanced. 
Dillon  and  Stoddard  (1977)  studied  the  heat  flux  to  turbine 
components  in  a  shock  tube.  Dillon  an  Nagamatsu  (1984)  studied  the 
rate  of  heat  transferred  to  the  walls  of  a  shock  tube  from  the  flow 
behind  a  shock  wave.  Novak  (1987)  investigated  heat  transfer  to  a 
flat  plate  by  using  thin  film  heat  flux  gages.  Jurgelwicz  (1989) 
developed  a  numerical  technique  to  obtain  heat  flux  from  the  data 
obtained  by  heat  flux  gages,  mounted  on  a  flat  plate,  exposed  to 
flow  behind  a  shock  wave.  He  also  extended  his  study  to  determine 
the  effects  of  film  cooling  on  the  rate  of  heat  transferred  to  a 
flat  plate  from  flow  behind  a  shock.  Rockwell  (1989)  developed  an 
analog  circuit  to  directly  record  the  heat  flux  measured  by  the 
thin  film  heat  flux  gages  in  terms  of  voltages.  This  direct 
conversion  of  voltages  into  heat  flux  tends  to  reduce  the  errors 
introduced  by  other  methods  like  numerical  techniques  etc. 

The  high  temperatures  encountered  in  modern  gas  turbines 
necessitated  studies  of  ways  of  protecting  exposed  surfaces  from  a 
hot  stream  of  gases.  One  such  way  is  to  inject  cooler  gas  under  the 
boundary  layer,  thus  forming  a  protective  film  along  the  surface. 
Goldstein  (1971),  in  his  book  "ADVANCES  IN  HEAT  TRANSFER",  has 
consolidated  the  work  of  many  researchers  pertaining  to  film¬ 
cooling.  It  provides  correlations  from  experimental  film-cooling 
studies  and  the  consequent  effects  of  film-cooling  on  the  rate  of 
heat  transfer. 

Inspite  of  the  commonalities  involved  between  film  cooling  in 
steady  and  transient  facilities  no  direct  comparison  between  the 
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two  can  be  made.  An  essential  element  to  understand  here  is  that 


the  flow  conditions  in  a  transient  facility  do  not  reach  a  steady 
state  at  any  time  and  are  thus  different  than  the  flow  conditions 
encountered  in  a  steady  state  facility  . 

Schlichting  (1987)  states  that  the  local  Nusselt  number  can  be 
increased  by  generating  free-stream  turbulence.  Blair  (1981)  found 
that  the  Reynolds  analogy  factor  increased  by  slightly  more  than 
one  percent  for  every  one  percent  increase  in  free-stream 
turbulence  intensity.  Simonich  and  Bradshaw  (1978)  determined  that 
grid  generated  free-stream  turbulence  increased  the  rate  of  heat 
transfer  by  about  five  percent  for  every  one  percent  increase  in 
the  longitudinal  turbulence  level. 

OBJECTIVES  AND  SCOPE 

This  Study  was  a  continuation  of  the  research  carried  out  by 
Jurgelwicz  (1989)  and  Rockwell  (1989).  the  general  purpose  of  the 
study  was  to  use  a  shock  tube  to  investigate  the  rate  of  heat 
transfer  to  a  flat  plate  under  the  transient  conditions  behind  the 
incident  shock  wave.  The  broad  objectives  of  this  study  were: 

1.  To  calibrate  the  heat  flux  gages  for  temperature 
coefficient  and  the  material  properties  forming  the  rooted  product 
of  density,  specific  heat  and  thermal  conductivity. 

2.  To  determine  the  rate  of  heat  transfer  to  a  flat  plate, 
with  round  leading  edge,  from  flow  behind  a  shock  wave  in  a  shock 
tube . 

3.  To  study  the  effects  of  free-stream  turbulence  on  the  rate 
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of  heat  transfer  to  a  flat  plate. 

4.  To  study  the  effects  of  film-cooling  on  the  rate  of  heat 
transfer  to  a  flat  plate. 

The  scope  of  the  research  ranged  from  generating  shock  Mach 
numbers  from  l.l  to  about  1.5  thereby  resulting  in  flow  behind  the 
incident  shock  wave  with  temperatures  and  pressures  as  high  as 
about  92  degrees  centigrade  and  60  in  of  Hg  respectively.  During 
the  film-cooling  analysis  the  blowing  ratio  of  the  film-cooling 
flow  was  varied  from  0.2  to  about  5.5  depending  on  the  flow 
conditions  of  the  primary  flow. 

An  understanding  of  the  phenomena  taking  place  during  the  heat 
transfer  to  a  flat  plate  and  the  effects  of  free-stream  turbulence 
and  film-cooling  on  the  rate  of  heat  transfer  would  help  in 
extending  the  acumen  to  the  complexities  involved  in  gas  turbine 
heat  transfer. 
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II .  THEORY 


THE  SHOCK  TUBE 

As  defined  by  Chapmann  And  Walker  (1971)  a  shock  tube  is  a 
device  that  allows  the  experimental  studies  of  many  of  the 
phenomena  associated  with  the  occurrence  of  shock  waves.  A 
conceptually  simple  shock  tube  is  shown  in  Figure  l.  It  essentially 
consists  of  two  sections  divided  by  a  diaphragm.  The  high  pressure 
section  of  the  tube  is  called  the  driver  section  whereas  the  low 
pressure  section  is  called  the  driven  section.  By  pressurizing  the 
driver  section  and  rupturing  the  diaphragm,  a  shock  wave,  the 
strength  of  which  is  determined  by  the  pressure  ratio  between  the 
two  sections,  is  generated.  This  results  in  a  discontinuity  in  the 
properties  of  the  medium  across  the  shock  wave.  This  can  best  be 
understood  with  the  help  of  Figures  2(a)  to  (d). 

At  time  t=0  the  high  pressure  gas  in  the  driver  section  is 
separated  from  the  low  pressure  gas  by  a  diaphragm.  When  the 
diaphragm  is  ruptured  a  shock  wave  travels  into  the  driven  section 
while  an  expansion  wave  travels  into  the  driver  section.  The 
initial  condition  in  the  driver  section  is  designated  as  region  4 
and  the  driven  section  as  region  1.  Before  the  two  waves  reach 
either  ends  of  the  shock  tube,  at  some  time  t-t^  the  region  notated 
as  region  3  in  Figure  2(a)  has  been  traversed  by  the  expansion  wave 
moving  into  the  driver  section  whereas  the  region  notated  as  region 
2  in  Figure  2(a)  has  been  traversed  by  the  shock  wave.  Though  the 
velocities  and  pressures  of  region  2  and  3  are  equal,  since  the  gas 
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state  in  region  3  is  obtained  by  an  isentropic  process  the  specific 
entropies  and  the  temperatures  of  the  two  regions  are  not  equal 
thus  implying  that  a  "contact  discontinuity"  exist  between  the  two 
regions.  This  contact  di scont inui ty  moves  in  the  same  direction  as 
the  shock  wave  with  a  velocity  and  pressure  V2=V3  and  P2=P3 
respectively  and  always  separates  the  shocked  region  2  from  the 
"expanded  region  3".  The  section  of  the  flow  with  uniform  pressure 
and  velocity,  i.e  region  2,  is  utilized  as  the  test  section  unless 
a  condition  with  high  temperature  and  high  pressure  at  zero 
velocity  is  desired  which  is  obtained  in  region  5  of  Figure  2(d). 

Figure  2(b)  shows  the  pressure  distribution  in  the  shock  tube 
at  t=tj  .  The  pressure  ratio  P4/P1  and  the  temperature  in  region  l 
determines  the  speed  of  the  shockwave  (Chapmann  and  Walker,  1971). 
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Note;  It  has  been  assumed  that  the  gas  in  the  driver  section 
is  the  same  as  that  in  the  driven  section. 

The  test  conditions  of  interest  for  this  investigation  are 
those  established  in  region  2.  The  properties  of  this  region  can  be 
determined  from  the  known  properties  of  region  1  and  the  strength 
of  the  shock  wave.  The  governing  equations,  extracted  from  the  text 
of  Shapiro  (1987),  to  determine  the  fluid  properties  behind  a 
normal  shock  wave  are: 
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Figure  2(c)  represents  the  temperature  distribution  withen  the 
shock  tube  at  time  t=tj  whereas  Figure  2(d)  essentially  indicates 
the  limited  time  available  for  a  test  to  be  carried  out  before  the 
properties  of  the  shocked  region  are  changed  by  a  reflected  shock 
wave  or  the  arrival  of  the  contact  discontinuity.  A  detailed 
account  of  the  shock  tube  phenomena  can  be  obtained  from  texts 
dealing  with  compressible  flow  dynamics. 


THE  BOUNDARY  LAYER. 

In  his  paper  "FLUID  MOTION  WITH  VERY  SMALL  FRICTION"  L. 
Prandtl  (1904)  proved  that  the  flow  about  a  solid  body  can  be 
divided  into  two  regions:  a  very  thin  layer  in  the  neighborhood  of 
the  body  where  friction  plays  an  important  role  and  a  region 
outside  this  layer  where  friction  may  be  neglected. 

This  very  thin  layer,  where  friction  plays  an  important  role, 
is  known  as  the  "boundary  layer".  The  problem  of  heat  flux  between 
a  solid  body  and  the  fluid  flow  past  it  belongs  to  the  class  of 
problems  where  boundary  layer  phenomena  plays  a  decisive  part 
( Schl i cht i ng ,  1987:  1-4). 
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Consider  a  round  leading  edge  flat  plate  suspended  in  a  shock 
tube.  When  an  incident  shock  wave  passes  over  the  flat  plate  a  thin 
boundary  layer  develops  over  it  as  a  result  of  the  gas  velocity 
established  in  region  2.  If  the  plate  temperature  is  T^,  the  shock 
velocity  is  Uj  and  the  temperature  and  velocity  of  the  free  stream 
flow  behind  the  shock  wave  is  T2  and  U2  then,  as  shown  in  Figure 
the  temperature  and  free  stream  velocity,  T2  and  U2  become  and 
zero  at  the  surface  of  the  plate. 

Blasius  (1908)  presented  solutions  to  equations  governing  the 
flow  in  the  steady  region  of  the  boundary  layer,  whereas  Mirels 
(1956)  obtained  empirical  solutions  to  equations  governing  the 
transition  region.  A  detailed  discussion  of  this  is  given  by  Davies 
and  Bernstein  (1969). 

HEAT  TRANSFER  THROUGH  THE  BOUNDARY  LAYER 

The  prime  focus  of  this  research  is  to  understand  the  process 
of  heat  transfer  and  the  influence  of  free-stream  turbulence  and 
film-cooling  on  heat  transfer  to  a  flat  plate.  The  pressure, 
temperature,  and  velocity  relationships  given  by  Equations  (3), 
(4),  (5),  used  in  conjunction  with  the  equations  presented  herewith 
provide  the  theoretical  heat  transfer  basis  for  comparison  with 
experimental  results. 

The  heat  transfer  to  the  flat  plate  in  a  high  speed  subsonic 
flow  is  given  by  (Kays  and  Crawford,  1980:299): 

^  ^aw~  (  6 ) 
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For  a  constant  property  perfect  gas,  the  adiabatic  wall 
temperature  is  given  by  the  relationship  (Mirels,  1956:23) 

of 

The  recovery  factor  "r"  for  a  steady  turbulent  boundary  layer 
can  be  approximated  as  (Mirels,  1956:23) 


The  heat  transfer  coefficient  h^^  can  be  determined  from  the 
expression  for  the  local  Nusselt  number  for  turbulent  flow  (Kays 
and  Crawford,  1980)  as: 


(9) 


This  would  then  imply  that  in  order  to  find  the  local  heat 
transfer  coefficient  we  need  to  know  the  local  Nusselt  number  which 
can  be  obtained  from  the  expression  (Kakac,  et  al ,  1987:  14-24) 


—  2 

NUj^=0 . 0 02 87i?ej/  Pr  ^  (10) 

The  local  Reynolds  number  can  be  obtained  from  the 
expression: 


i?e=3f 


(11) 


The  Prandtl  and  Stanton  numbers  for  a  steady  turbulent  flow 
are  expressed  as: 


5tPr '’  •* =0 . 02  87 •  2 
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(13) 


Since  the  fluid  properties  are  not  constant  across  the 
boundary  layer  in  high  speed  flow,  they  need  to  be  evaluated  at 
some  reference  temperature  T* .  According  to  Eckert  (1976),  the  low 
speed  constant  properties  correlations  for  Nusselt  number  can  be 
used  for  air  for  Mach  numbers  up  to  20  with  an  approximate  error 
not  exceeding  more  than  just  a  few  percent. 

The  expression  for  reference  temperature,  as  given  by  Kays  and 
Crawford  (1980)  is: 


T*  =  Is!^^0 . 22 


Equations  (6)  to  (14)  constitute  the  theoretical  results  for 
the  steady  turbulent  boundary  layer. 


ELECTRICAL  ANALOG  FOR  HEAT  TRANSFER 


One  dimensional  heat  transfer  to  a  semi  infinite  medium  is 
expressed  as  : 

dx  ^ 

(16) 

Equations  (15)  and  (16)  can  be  reduced  to  a  diffusion  equation 
by  substitution  to  get: 


dT_  k 
dt  pCp  dx^ 
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The  solution  to  Equations  (15),  (16)  and  (17)  for  a  semi¬ 
infinite  medium  is  given,  in  the  Laplace  domain,  by  OLdfield  (1978) 
as : 


g=^p(^jjc.y/s.  X  (13) 

One  dimensional  diffusion  of  electrical  charge  through  a 
medium  with  distributed  capacitance  per  unit  volume  and  resistance 
per  unit  length  is  given  by  Oldfield  (1978)  as: 


dx  dt 


(1^) 


1  dV 
r  dx 


Combining  Equations  (19)  and  (20)  gives: 


dv_  1 

dt  TC 


(n) 


Equations  (19),  (20),  and  (21)  can  be  solved  in  the  Laplace 


domain  to  get: 


Note  that  Equations  (18)  and  (22)  have  identical  forms.  As 
such  analogous  thermal  and  electrical  systems  can  be  constructed. 
Theoretical  insight  for  analogy  and  construction  of  the  analogous 
electrical  system  for  heat  flux  is  documented  by  Roclcwell  (1989). 
Only  the  final  form  of  the  heat  flux  equation,  as  obtained  by 
Oldfield  (1984)  is  presented  here  as: 
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(2.) 

The  term  V ( pCpk)  is  a  calibration  constant  for  the  heat  flux 
gage,  V^ak  is  the  temperature  coefficient  of  the  heat  flux  gage  and 
( 1/G)^ ( r/c )  is  the  calibration  constant  for  the  heat  flux  analog 
circuits.  The  calibration  constants  for  the  heat  flux  analog 
circuits  used  in  this  study  were  obtained  by  the  Rockwell  and  are 
reproduced  as  Table  I  .  For  details  please  refer  to  Rockwell 
( 1989)  . 


Table  I:  Calibration  Data  for  the  Heat  Flux  Analog  Circuit. 


CIRCUIT 

SLOPE 

ST.  DEVIATION 

CORRELATION 

31-105 

0.74371 

1.928 

0.99945  1 

31-790 

0.78699 

1.957 

0.99942  1 

31-820 

0.73212 

1.670 

0.99958  1 

31-850 

0.73858 

2.000 

0.99941  1 

31-870 

0.79757 

1.248 

0.99955 

32-100 

0.73695 

2.018 

0.99939 

32-200 

0.70370 

2.149 

0.99927  1 

HEAT  TRANSFER  WITH  FREE -STREAM  TURBULENCE 

Free-stream  turbulence  is  defined  as  the  ratio  of  velocity 
fluctuations  about  the  mean  velocity  to  the  mean  flow  velocity: 
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(24: 


Free-stream  turbulence  serves  to  increase  the  rate  of  heat 
transfer.  This  increase  in  the  rate  of  heat  transfer  results  from 
the  fact  that  the  free-stream  eddies  get  entrained  in  the  boundary 
layer  thereby  transitioning  the  flow  from  laminar  to  turbulent.  The 
free-stream  eddies  also  have  a  profound  effect  in  increasing  the 
heat  transfer  in  the  case  where  the  boundary  layer  is  already 
turbulent . 

Simonich  and  Bradshaw  (1978)  found  that  the  ratio  of  Stanton 
numbers  of  flows  with  and  without  turbulence  can  be  established  for 
low  turbulence  levels  (Simonich  and  Bradshaw,  1978  :67l) 


St 


=l+5ru 


(25) 


HEAT  TRANSFER  WITH  FILM-COOLING 

As  Stated  earlier,  film  cooling  is  an  effective  way  of  forming 
a  protective  film  on  a  surface  by  injecting  a  cooler  fluid  under 
the  boundary  layer  which  is  established  by  the  flow  of  a  hot  stream 
of  air  over  the  surface. 

The  injected  fluid  can  enter  the  boundary  layer  in  a  number  of 
different  ways.  However,  to  keep  the  analysis  relatively  simple  a 
flat  plate  with  injection  holes  at  90°  to  the  main  hot  stream  of 
air  has  been  used  during  this  study. 

With  the  assumption  that  film  cooling  produces  turbulent 
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boundary  layer  (Kakac,  et  al,  1987:2-28),  the  heat  transfer  rate  Is 
expressed  as: 

^  ^aw~  (  2  0 ) 

It  should  be  kept  in  mind  that  the  adiabatic  wall  temperature 
for  flow  with  film  cooling  is  different  than  the  adiabatic  wall 
temperature  for  flow  without  film  cooling. 

Though  a  number  of  studies  on  film  cooling  has  been  done  in 
steady  state  facilities  like  wind  tunnels,  etc;  no  such  literature 
is  available  for  studies  carried  in  transient  facilities  like  shock 
tube.  The  most  intriguing  problem  in  a  transient  facility  like  the 
shock  tube  is  that  the  flow  does  not  get  sufficient  time  to  attain 
steady  state  conditions  (steady  wall  temperature).  This  results  in 
a  situation  where  the  adiabatic  wall  temperature  cannot  be  achieved 
and  thus  measured. 

The  experimental  data  for  this  part  of  the  study  can  thus  only 
be  presented  in  terms  of  physically  measurable  quantities  like  the 
heat  flux,  blowing  rate,  and  the  Reynolds  number  without  any 
theoretical  or  experimental  comparison. 
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Ill  EXPERIMENTAL  APPARATUS 


SHOCK  TUBE 

The  shock  tube  used  for  this  study  was  the  AFIT  low  pressure 
shock  tube.  This  shock  tube  has  four  main  sections:  the  driver 
section,  the  driven  (expansion)  section,  the  10.16  cm  x  20.32  cm 
test  section  and  the  dump  tank  section.  Figure  4  shows  a  schematic 
of  the  shock  tube.  The  high  and  low  pressure  sections:  that  is,  the 
driver  and  the  driven  sections  are  separated  by  a  Mylar  diaphragm 
(available  in  thicknesses  varying  from  25.4  um  to  177.8  ym)  their 
use  being  dependent  on  the  required  strength  of  the  shock  wave. 
For  most  of  this  study  127.0  jim  thick  diaphragms  were  used. 
Ordinary  compressed  dry  air  of  up  to  100  psig  was  used  to 
pressurize  the  driver  section  of  the  shock  tube.  The  pressure  of 
the  driver  section  "P^"  was  measured  by  a  calibrated  pressure  gage 
in  inches  of  mercury.  The  pressure  of  the  driven  section  "Pj"  was 
assumed  to  be  the  same  as  ambient  pressure.  The  diaphragm  was 
ruptured  by  a  mechanically  controlled  plunger. 

INSTRUMENTED  FLAT  PLATE 

An  instrumented  flat  plate  was  used  as  the  test  model  for  the 
heat  transfer  study  and  was  installed  in  the  center  of  the  test 
section.  The  plate  was  64.008  cm  long,  10.16  cm  wide,  1.915  cm 
thick  and  had  a  round  leading  edge.  It  was  instrumented  with  seven 
heat  flux  gages  and  a  single  row  of  forty  one  film-cooling  holes. 
Each  film  cooling  hole  was  1  mm  in  diameter  and  separated  1.9812  mm 
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center  to  center.  Figure  5  illustrates  the  instrumentation  of  the 
flat  plate. 

FILM  COOLING  SYSTEM 

The  film-cooling  holes  were  located  5.08  cm  downstream  of  the 
leading  edge  of  the  flat  plate.  These  holes  were  supplied  with  dry 
filtered  air  from  the  main  air  supply  system.  The  air  pressure  was 
regulated  through  a  regulating  system  consisting  of  a  dome  valve 
and  a  pressure  regulated  air  bottle.  Bottle  air  of  known  pressure 
was  applied  to  the  reference  side  of  the  dome  valve  through  an 
electrically  controlled  solenoid  valve.  This,  in  turn,  controlled 
the  pressure  of  the  dry  filtered  air  (approximately  the  same 
pressure)  to  a  chamber  supplying  air  to  the  film-cooling  holes.  The 
air  pressure  for  the  row  of  film-cooling  holes  was  measured  by  a 
pressure  transducer  located  in  the  pipe  connected  to  the  air 
chamber.  Figure  6  shows  a  schematic  of  the  film-cooling  system. 

TURBULENCE  GENERATING  SYSTEM 

A  turbulence  generator  system  was  used  to  form  an  air  grid 
across  the  shock  tube  so  as  to  generate  free-stream  turbulence  in 
the  flow  over  the  flat  plate.  The  turbulence  generator  was 
installed  just  in  front  of  the  leading  edge  of  the  flat  plate. 
Figure  6  illustrates  a  schematic  of  the  turbulence  generating 
system.  Regulated,  dry,  filtered  air  was  injected  perpendicular  to 
the  flow  at  twelve  locations:  two  each  on  top  and  bottom  walls  of 
the  turbulence  generator  and  four  each  on  the  side  walls.  The 
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supply  of  high  pressure  air  at  these  locations  formed  a  strong  air 
grid  across  the  on-coming  flow  in  the  shock  tube.  This  yielded  an 
increase  in  the  free-stream  turbulence  intensity. 

FLOW  VISUALIZATION 

The  schlieren  flow  visualization  system  was  used  to  obtain 
images  of  the  boundary  layer  and  shock  patterns  across  the  test 
section.  Two  mirrors  of  30"  focal  length,  along  with  a  spark 
source,  a  camera  and  a  knife  edge,  were  used  to  obtain  photographs 
of  the  flow  across  the  plate.  The  spark  source,  supplied  by  a  high 
voltage  power  supply,  was  triggered  by  a  delayed  voltage  signal 
from  the  front  pressure  transducer.  The  voltage  signal  was  delayed 
by  a  combination  of  the  Cordin  proportional  delay  generator  Model 
435  and  a  Cordin  delay  generator  Model  453.  The  flow  images  were 
captured  on  Polaroid  type  52  (ASA  3000)  film.  Figure  7  illustrates 
the  configuration  of  the  schlieren  system. 

INSTRUMENTATION 

1.  THIN- FILM  HEAT  FLUX  GAGES 

Seven  thin-film  heat  flux  gages  were  mounted  along  the 
centerline  of  the  plate,  spaced  proportionally  downstream  of  the 
film  cooling  holes.  Figure  8  shows  the  proportional  location  of  the 
heat  flux  gages  along  the  instrumented  plate.  The  thin-film  gages 
are  made  of  a  platinum  film  0.5  mm  wide  and  0.1  ym  thick  deposited 
on  a  Pyrex  7740  substrate.  The  /  (pCpk)  value  of  the  heat  flux 
gages  was  obtained  by  calibration  of  the  gages.  The  temperature 
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coefficient  "V^ak"  used  in  equation  3  was  also  obtained  by 
calibration.  The  heat  flux  gages  available  for  this  study  were  gage 
numbers  2,  4,  and  6,  counting  from  the  leading  to  trailing  edge  of 
the  plate.  These  gages  were  powered  by  an  excitation  voltage  from 
the  PSC  8015-1  dc  amplifier. 

2.  PRESSURE  TRANSDUCERS 

Two  Endevco  Model  8530A-100  pressure  transducers  were 
flush  mounted  on  the  top  wall  of  the  shock  tube.  The  forward 
pressure  transducer  was  1.1811  m  upstream  of  the  leading  edge  of 
the  flat  plate  whereas  the  rear  pressure  transducer  was  located 
0.4699  m  upstream  of  the  plate  leading  edge.  A  third  pressure 
transducer  was  located  upstream  of  the  film  cooling  air  supply  line 
to  ensure  accurate  measurements  of  the  film  cooling  air  pressure. 
The  excitation  voltage  to  all  three  pressure  transducers  was 
supplied  by  a  Power  Mate  Corporation  power  supply.  The  main  purpose 
of  the  pressure  transducers  in  the  shock  tube  was  to  measure  the 
time  interval  between  their  responses  to  the  pressure  change 
induced  by  the  passage  of  a  shock  wave.  This  lapse  time  was  then 
used  to  calculate  the  speed  of  shock  wave  and  thus  its  Mach  number. 
During  schlieren  photography  the  output  of  the  front  and  rear 
pressure  transducers  was  supplied  to  the  time  delay  generators  to 
control  the  timing  for  the  photography.  The  front  pressure 
transducer  was  also  used  to  trigger  the  DL-1200  Datalab  recorder, 
to  be  described  forthwith.  Calibration  details  of  these  pressure 
transducers  is  given  in  Appendix  D. 
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3.  HOT  WIRE  ANEMOMETER 


To  measure  the  turbulence  level  of  the  flow  behind  the 
normal  shock  wave,  a  TSI  model  intelligent  flow  analyzer  model  IFA 
100,  embedded  with  a  hot  wire  anemometer,  was  used  to  provide  and 
measure  the  hot  wire  bridge  output  voltage.  A  TSI  model  1214-10  hot 
film  sensor  was  mounted  in  the  test  section  to  collect  the  required 
data . 


4.  THERMOCOUPLE 

The  calibration  procedure  to  obtain  the  temperature 
coefficient  of  the  heat  flux  gages  required  a  thermocouple  to 
measure  the  temperature  of  the  water  bath.  A  K  type  (chromel- 
alumel)  thin  foil  thermocouple  was  used  in  conjunction  with  a 
battery  powered  ice  point  for  temperature  measurements.  The 
temperature  of  the  water  bath  was  indicated  in  °F  on  a  previously 
calibrated  Omega  Digicator. 

5.  SIGNAL  GENERATOR 

The  procedure  used  to  obtain  the  V  (pCpk)  values  of  the 
heat  flux  gages  required  that  a  square  pulse  of  known  amplitude, 
period  and  width  be  applied  to  the  heat  flux  gage  connected  as  one 
leg  of  a  bridge  circuit.  Application  of  such  a  pulse  raised  the 
temperature  of  the  gage.  This  changed  the  resistance  of  the 
platinum  foil  in  the  heat  flux  gage  which  in  turn  produced  a 
measurable  voltage  by  the  bridge  circuitry.  To  produce  a  single 


21 


square  pulse  of  known  amplitude,  period  and  width  a  Wavetek  Model 
signal  generator  was  used  during  the  calibration  phase. 

6.  PSC  MODEL  BRIDGE/AMPLIFIER 

The  response  of  a  heat  flux  gage  to  any  temperature 
change  can  be  measured  by  replacing  one  leg  of  a  Wheatstone  bridge 
circuit  with  the  heat  flux  gage.  A  temperature  change  in  the 
environments  of  the  heat  flux  gage  unbalance  the  previously 
balanced  bridge  and  thus  the  output  of  the  bridge.  The  PSC  8115 
provided  the  three  legs  of  the  bridge  circuitry  with  precision 
resistors  and  a  variable  resistor  in  parallel  to  the  leg  opposite 
to  the  heat  flux  gage,  thereby  enabling  the  user  to  balance  the 
bridge.  The  constant  bridge  excitation  voltage  was  set  at  2.5  volts 
for  each  bridge  circuit  in  series  with  a  1000  n  precision  resistor. 
This  limited  the  supply  current  to  a  nominal  value  of  2.5  mamps  to 
avoid  damage  to  the  heat  flux  gages. 

Any  temperature  change  sensed  by  the  heat  flux  gage  changed 
the  resistance  of  the  gage  by  a  proportional  amount  and  thus 
unbalanced  the  bridge.  The  output  of  the  bridge,  being  extremely 
low,  required  high  amplification.  Thus  an  amplifier  circuit  was 
also  required  for  each  bridge  module  in  use.  PSC  8015-1  high  gain 
differential  amplifiers  were  used  for  this  purpose.  The 
bridge/amplifier  circuits  for  each  heat  flux  gage  were  card  mounted 
and  installed  in  a  Transamerica  instruments  rack  designed  for  these 

modules . 
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6.  HEAT  TRANSFER  ANALOG  CIRCUIT 


The  heat  transfer  analog  circuit,  designed  by  Rockwell 
(1989)  was  used  to  convert  the  amplified  output  of  the  heat  flux 
gages  into  potential  heat  flux.  The  output  of  the  heat  flux  analog 
in  combination  with  the  calibration  constants  resulted  in  the 
measured  heat  flux  (Equation  31).  A  block  diagram  of  the  analog 
circuit  is  shown  in  Figure  9.  The  details  of  construction  and 
working  of  the  heat  flux  analog  circuit  were  reported  by  Rockwell 
( Rockwel  1 ,  1989 )  . 

7.  WAVE  FORM  RECORDER 

A  high-speed  digital  data  recorder,  Datalab  DL  1200,  was 
used  to  record  the  required  output  from  the  respective  instruments. 
The  DL-1200  waveform  recorder  is  a  12  bit  device  that  transforms  an 
analog  voltage  input  to  a  digital  value  for  eight  channels 
simultaneously.  Each  channel  has  its  own  gain  and  filter  settings. 
During  the  course  of  this  study  each  channel  was  sampled  at  its 
maximum  sampling  rate,  i.e  one  sample  per  channel  every  2  us.  The 
recorded  data  from  each  test  was  then  downloaded  from  the  DL-1200 
12  bit  memory  to  Zenith  16  bit  computer  and  stored  for  analysis. 

8.  SIGNAL  DELAY  GENERATORS 

The  Cordin  proportional  and  simple  delay  generators  were 
used  to  delay  the  voltage  signals,  from  the  front  pressure 
transducer,  used  to  trigger  the  spark  source  of  the  schlieren 
system.  The  delay  time  on  both  the  delay  generators  can  be  set 
such  that  the  combination  of  the  two  equals  the  time  required  for 
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the  shock  wave  to  be  positioned  at  a  point  desired  by  the  operator. 
The  combination  of  two  delay  generators  was  relatively  easier  to 
handle  than  using  Cordin  Model  453  delay  generator  only. 
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IV  EXPERIMENTAL  PROCEDURE 


The  basic  requirement  for  any  experimental  set-up  is  to  have 
a  sensing  element  which  can  discern  the  requisite  information  for 
further  analysis.  As  mentioned  earlier,  the  sensing  elements  used 
during  this  study  to  measure  heat  flux  were  the  thin  film  heat  flux 
gages.  However,  to  analyze  the  data  obtained  by  the  heat  flux 
gages,  knowledge  about  some  of  the  properties  of  the  heat  flux 
gages  was  required.  These  properties  were  the  temperature 
coefficient,  V^ak,  and  bulk  thermal  diffusivity,  V  pCpk.  These 
properties  can  either  be  obtained  from  theoretical  data  or 
established  by  calibration.  To  acquire  credible  results  it  was 
desirable  to  choose  the  latter  option  of  calibrating  the  heat  flux 
gages . 

CALIBRATION  OF  HEAT  FLUX  GAGES 

1.  CALIBRATION  FOR  TEMPERATURE  COEFFICIENT 

The  calibration  circuit  used  to  obtain  the  temperature 
coefficient  of  the  heat  flux  gages  is  shown  in  Fig  10.  The 
calibration  procedure  adopted  was  as  follows:- 

The  heat  flux  gages,  flush  mounted  on  the  flat  plate,  were 
loosened  such  that  with  the  plate  inverted  they  could  dangle 
freely. 

The  gage  to-be-calibrated  was  then  connected  as  one  leg  of  the 
Wheatstone  bridge  already  built  in  the  PSC  8115  bridge  module. 
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The  PSC  8015-1  differential  amplifier  controls  were  set  at  a 
gain  250  and  a  filter  setting  of  10  kHz.  The  circuit  was  then 
zeroed  out  by  cycling  through  the  internal  zeroing  routine  pre¬ 
programmed  into  the  module. 

Setting  the  bridge/amplifier  module  at  an  AMP-D  mode  and 
connecting  the  bridge  output  to  a  previously  zeroed  voltmeter,  the 
variable  resistor  on  the  PSC  8115  bridge  module  was  used  to  balance 
the  bridge. 

The  heat  flux  gage  was  then  placed  in  a  water  proof,  thin 
plastic  bag  along  with  a  J-type  thermocouple.  An  Omega  digicator, 
which  sensed  the  thermocouple  output,  referenced  it  to  its 
internal  ice  point  and  displayed  the  temperature  of  the  water  bath 
in  degrees  Fahrenheit. 

A  1000  ml  beaker  filled  with  about  750  ml  of  water  was  placed 
on  a  thermal  mixer.  A  magnetic  stirrer  was  then  placed  inside  the 
beaker.  The  plastic  bag,  containing  the  heat  flux  gage  and  the  J- 
type  thermocouple,  was  then  submersed  in  the  water  bath  such  that 
the  water  could  not  enter  the  bag. 

The  thermal  mixer  was  turned  "on"  so  as  to  slowly  mix  the 
water  to  eliminate  any  hot  spots.  The  heating  element  was  then 


turned  "on" 

to  a  slow 

heating  position 

so 

as  to  minimize 

the 

di f  f erence 

between  the 

response  times 

of 

the  gage  and 

the 

thermocouple.  A  schematic  of  the  instrument  interfacing  for  this 
calibration  is  shown  in  Figure  li. 

A  change  in  the  digicator  display  was  noted  along  with  the 
output  voltage  of  the  now  unbalanced  bridge.  A  total  of  li  such 
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readings  of  the  digicator  and  the  output  voltage  of  the  bridge, 
indicating  the  response  of  the  heat  flux  gage  to  the  change  of 
temperature  of  its  environments,  were  taken  for  each  heat  flux 
gage. 

The  calibration  data  was  then  plotted  and  the  slope  of  the 
fitted  curve  gave  the  temperature  coefficient  of  the  heat  flux 
gage.  Table  II  shows  the  temperature  coefficients  to  the  heat  flux 
gages  calibrated  for  this  study.  The  calibration  plots  are  shown  in 
Appendix  C. 


Table  II:  Temperature  Coefficients  of  the  Heat  Flux  Gages. 


Gage  Number 

Gage  Serial  Number 

a  (V/°F) 

2 

516 

0.0241006 

4 

503 

0.0236467  I 

6 

120 

0.0178865  1 

1 

2.  CALIBRATION  FOR  BULK  THERMAL  DIFFUSIVITY  V  (oCpK) 

Thin  film  heat  flux  gages  consist  of  a  very  thin  platinum 
film,  0.5  mm  wide  and  0.1  ym  thick,  deposited  on  Pyrex  7740 
substrate.  The  calibration  method  required  to  obtain  the  V  (pCpk) 
values  of  the  Pyrex  7740  substrate  used  in  the  heat  flux  gages  is 
a  development  of  a  pulse  technique.  Figures  12  and  13  show  the 
bridge  circuit  and  the  instrument  interfacing  used  for  calibration. 

A  step  current,  of  known  amplitude,  is  applied  to  the  heat 
flux  gage  which  acts  as  one  leg  of  the  balanced  bridge.  The 
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passage  of  current  allows  the  platinum  film  to  heat  up  thereby 
changing  its  resistance.  This  change  in  resistance  unbalance  the 
bridge  and  the  output  voltage  can  be  recorded  to  the  DL-1200. 

A  small  droplet  of  liquid,  of  known  bulk  thermal  diffusivity, 
(in  this  case  glycerine)  is  then  deposited  on  the  surface  of  the 
heat  flux  gage.  With  the  bridge  balanced,  current  of  the  same 
amplitude  as  before  is  again  passed  through  the  heat  flux  gage. 
This  time  heat  flows  both  into  the  substrate  of  the  gage  and 
glycerine.  The  output  voltage  of  the  now  unbalanced  bridge  is 
recorded  to  the  DL-1200. 

The  two  voltage  records  have  the  same  time  dependence  but 
different  amplitudes.  For  illustration  refer  to  Figure  14.  As  can 
be  seen  from  the  figure,  the  output  voltage  of  both  the  experiments 
is  parabolic  in  nature.  This  would  imply  that  a  linear  relationship 
can  be  established  between  the  change  in  resistance  aR  of  the 
platinum  f'^11  of  heat  flux  gages  and  iTT. 

r  ■ 

If  aRj/1  t  is  the  slope  of  the  linearly  fitted  curve  obtained 
from  an  experiment  in  air  and  AR2//~t  is  the  slope  of  the  linearly 
fitted  curve  obtained  from  an  experiment  in  glycerine  then, 
following  the  analysis  given  in  Shultz  and  Jones  (1973),  the 
correct  value  of  J  pCpk  of  the  Heat  Flux  gage  substrate  can  be 
obtained  from  the  relationship: 


y/t 

^1/ 


{-) 


A  complete  derivation  of  this  expression  is  given  as  Appendix 
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A. 

Since  constant  current  of  equal  amplitudes  is  passed  through 
the  heat  flux  gages  in  both  experiments.  Equation  (27)  can  also  be 
written  as: 


^^^=9 25 sec  *  i2  j) 

Table  III  shows  the  results  obtained  during  the  calibration  of  the 
thin-film  heat  flux  gages  used  to  carry  out  this  study. 


Table  III:  Bulk  Thermal  Diffusivity  of  Heat  Flux  Gages. 


Gage  Number 

516 

0. 101388 

0.063046 

1521 

503 

0.081732 

0.052482 

1660 

120 

0.017039 

0.010646 

1540 

A  detailed  calibration  procedure  for  determining  the  bulk 
thermal  diffusivity  of  the  heat  flux  gages  used  during  this 
research  is  given  in  Appendix  B.  The  calibration  curves  are  also 
given  in  Appendix  B. 
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SHOCK  GENERATION 


To  achieve  repeatable  free-stream  flow  conditions  a  Mylar 
diaphragm  of  the  same  thickness  was  used  and  the  shock  speed  was 
monitored  for  every  run.  The  same  shock  speed  for  two  different 
runs  essentially  means  that  the  free-stream  flow  properties  behind 
the  shock  wave  are  same  for  both  the  runs.  This  implies  then  that 
they  would  also  produce  the  same  heat  transfer  at  some  given  point 
on  the  flat  plate.  During  the  film-cooling  and  turbulence 
generation  runs,  to  avoid  any  flow  anomalies  from  occurring  in  the 
boundary  layer,  film-cooling/turbulence  generation  were  initiated 
at  approximately  the  same  time  the  shock  was  generated  by  bursting 
the  Mylar  diaphragm  between  the  high  and  low  pressure  sections  of 
the  shock  tube. 

DATA  COLLECTION 

All  the  data  required  for  analysis  was  amplified  and  filtered 
through  the  PSC  81015-1  differential  amplifier.  The  amplifier  was 
set  at  a  gain  of  1000  and  a  Low  Pass  filter  of  3  kHz.  This 
amplified  and  filtered  data  was  then  converted  into  heat  flux 
voltages  through  the  heat  flux  analog  and  recorded  in  the  datalab 
DL-1200  recorder.  The  DL-1200  recorder  was  remotely  controlled  by 
an  upgraded  Z-248  computer.  For  almost  all  the  test  runs  the  DL- 
1200  was  set  at  the  following  settings: 


B - 

1  TRIGGER  LEVEL 

PRE  DELAY 

SAMPLE  INTERVAL  I 

fl 

8 

! 

6% 

2  /isec  i 

■ 
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A  NEWDL  software  package  was  used  to  interface  the  DL-1200  and 


Z-248  computer  whereas  a  POST  program,  developed  by 
(1989),  was  used  for  data  reduction  and  analysis. 


Rockwel 1 
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V.  DATA  REDUCTION 


SHOCK  MACH  NUMBER 

The  shock  Mach  number  was  measured  for  every  run.  The  passage 
of  a  shock  across  the  front  and  rear  pressure  transducers  and  the 
heat  flux  gages  caused  a  significant  increase  in  the  output  voltage 
of  each  sensor  because  of  changed  flow  conditions.  The  time  of 
passage  of  the  shock  wave  across  these  sensors  can  be  determined 
from  their  output.  Knowing  the  lapse  time  between  the  passage  of 
the  shock  wave  across  the  front  and  rear  pressure  transducers  and 
the  distance  between  them,  the  shock  speed  was  evaluated  by  using 
the  equation: 


(30 


The  temperature  Tj  of  the  driver  section  of  the  shock  tube  is 
assumed  to  be  the  same  as  ambient  temperature.  The  Mach  number  of 
the  shock  wave  was  determined  by  using  equation: 


(31) 


Having  calculated  the  shock  Mach  number,  the  temperature, 
pressure,  Mach  number, and  the  density  of  the  flow  behind  the  shock 
wave  and  its  velocity  with  respect  to  the  hardware  of  the  shock 
tube  was  determined. 

TURBULENCE  INTENSITY 

The  purpose  of  measuring  the  free-stream  turbulence  and  the 
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generated  turbulence  intensities  was  to  verify  the  measurements 
made  by  Rockwell  (1989)  and  study  their  effect  on  the  heat  transfer 
solution  of  the  flat  plate. 

The  procedure  for  turbulence  measurements  was  also  adopted 
from  Rockwell  (1989).  Several  test  runs  were  made  using  the  same 
driver  pressure  of  100  in  of  Hg  but  different  hot  film  overheat 
parameter.  The  TSI  1214-10  hot  film  probe  with  a  90°  arm  and  probe 
support  was  placed,  through  the  top  wall  of  the  shock  tube,  such 
that  the  hot  film  was  perpendicular  to  the  flow  and  about  3  cm 
above  the  forth  heat  flux  gage.  The  90®  arm  support  was  used  to 
minimize  the  vibrations  caused  by  the  passage  of  the  shock  wave. 
Beside  the  TSI  1214-10  hot  film  probe  other  probes  like  TSI  1214- 
T1.5  hot  wire  probe  and  TSI  1214-20  hot  film  probes  were  also 
tried.  However,  as  can  be  seen  from  Figure  15,  the  TSI  1214-T1.5 
hot  wire  probe  showed  periodic  oscillations  caused  by  the 
vibrations  due  to  the  passage  of  the  shock  wave  whereas  the  TSI 
1214-20  hot  film  probe  was  too  slow  to  respond  during  the  short 
duration  of  the  test  runs.  High  frequency  vibrations  could  also  be 
seen  in  the  data  obtained  from  the  TSI  1214-10  hot  film  probe  but 
since  most  of  the  turbulence  data  occurs  between  5  Hz  and  10  kHz 
frequencies  above  10  kHz  were  filtered  out  without  damaging  the 
integrity  of  the  data. 

The  TSI  1214-10  hot  film  probe  was  connected  to  the  IFA  100. 
The  output  of  the  IFA  100  was  recorded  on  DL  1200  which  was 
triggered  by  a  signal  from  the  front  pressure  transducer  generated 
by  the  passage  of  a  shock  wave. 
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The  POST  program  was  used  to  filter  out  the  high  frequency 
spikes  from  the  data  and  then  used  to  calculate  the  mean  and  rms 
values  of  the  data.  The  turbulence  sub-routine  in  the  POST  program 
calculated  and  plotted  a  least  quadratic  fit  between  the  ratio 
^ film  over  heat  parameter.  The  turbulence 
intensity  of  the  flow  was  obtained  from  the  y-intercept  of  the 
least  quadratic  curve  fit.  For  theoretical  details  refer  to 
Rockwel  1  ( 1989  )  . 

The  same  procedure  was  repeated  with  the  turbulence  generator 
"on"  and  the  results  were  then  compared  with  those  obtained  by 
Rockwel  1  ( 1989  )  . 

HEAT  TRANSFER 

Figure  16  shows  the  circuitry  and  interfacing  of  the 
instrumentation  used  to  determine  the  heat  flux  into  the  flat 
plate. 

The  flow  behind  the  shock  is  at  an  elevated  temperature 
compared  to  the  flat  plate  and  the  flow  in  front  of  the  shock  wave. 
This  flow  of  hot  air  over  the  plate  results  in  heat  transfer  to 
the  flat  plate.  The  heat  flux  gages  sense  the  change  of  the  flat 
plate  temperature  caused  by  the  change  in  its  environment  and 
respond  almost  instantaneously.  The  change  in  resistance  of  the 
platinum  foil  of  the  heat  flux  gages  unbalance  their  respective 
bridge  circuitry  in  the  PSC  8115  modules.  This  results  in  a  voltage 
output  by  the  PSC  8115  modules  which  is  amplified  and  filtered  by 
the  differential  amplifiers  in  the  PSC  8105-1  modules.  The 
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ampl i f i ed  and  filtered  output  of  the  PSC  8015-1  modules  is  fed  into 
the  heat  flux  analog  which  converts  the  raw  voltages  into  heat  flux 
voltages.  These  voltages  are  recorded  by  the  DL-1200  for  storage  by 
the  Z-248  computer.  Multiplying  the  heat  flux  voltages  with  the 
calibration  constants  of  their  respective  heat  flux  gages  give  the 
actual  heat  transfer  into  the  flat  plate  at  each  specific  location. 


FILM-COOLING 

The  process  of  data  collection  and  reduction  for  film-cooling 
runs  was  essentially  the  same  as  that  for  the  heat  transfer  runs. 
The  main  feature  of  the  film-cooling  runs  was  to  determine  the 
blowing  ratio,  B,  and  its  consequent  effects  on  the  rate  of  heat 
transfer.  The  pressure  of  the  coolant  flow  was  determined  from  the 
pressure  transducer  located  a  short  distance  from  the  coolant 
chamber  of  the  flat  plate.  Refer  to  Figure  6  for  an  illustrative 
schematic  of  the  film  cooling  system. 

Knowing  both  the  pressures  of  the  coolant  flow  and  the  flow 
behind  the  shock  wave  the  velocity  of  the  coolant  flow  can  be 
determined  from  equation: 


=  +  (3v) 

Since  the  velocity  and  the  density  of  the  coolant  and  the 
primary  flows  are  now  known,  the  blowing  ratio,  B,  into  the  primary 
flow  can  be  determined  from  the  relationship  expressed  as: 
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(33) 


B=' 


PqVc 

P2V2 


If  the  pressure  of  the  coolant  flow  is  continuously  increased 
a  stage  would  come  when  the  Mach  number  of  the  coolant  flow  would 
reach  unity.  The  velocity  of  the  flow  at  this  point  can  be 
determined  by  the  equation  ( Zucrow  and  HOffman,  1976): 


2  v-5 
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Under  such  conditions  the  flow  is  said  to  be  choked.  Any  further 
increase  in  the  pressure  of  the  coolant  flow  would  not  increase  the 
velocity  of  the  flow.  Instead,  the  density  of  the  flow  would 
continue  to  increase  and  thus  increase  the  blowing  rate.  The 
density  of  the  choked  coolant  flow  can  be  determined  from  the 
equation : 


Pc“^oc' 
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VI.  RESULTS  AND  DISCUSSION 


SHOCK  MACH  NUMBER 

The  theoretical  Mach  number  of  a  shock  wave  can  be  calculated 
by  using  Equation  2.  However,  experimental  observations  show  that 
for  identical  pressure  ratio  P4/P1  the  shock  Mach  number  calculated 
experimentally,  using  Equation  (31),  is  not  quite  the  same  as  that 
calculated  theoretically.  Table  IV  shows  the  minor  variations 
between  the  experimental  and  theoretical  shock  Mach  numbers  for  the 
same  pressure  ratio. 


TABLE  IV:  COMPARISON  OF  MEASURED  AND  THEORETICAL  SHOCK  Mach  Number 


DATA  RUN 

NUMBER 

PRESSURE 

P4”  Hg 

PRESSURE 

Pf  Hg 

EXPERIMENTAL 

SHOCK  MACH  NO 

THEORETICAL  1 

SHOCK  MACH  NO  | 

R002 

60 

29.07 

1.209 

1.268  1 

R004 

80 

29,07 

1.269 

1.323  1 

R006 

100 

29.07 

1.320 

1.369  1 

R008 

120 

29.07 

1.360 

1.409  1 

I  FC603 

60 

29.30 

1 .  197 

1.260  1 

1  FC605 

60 

29.30 

1 .  197 

1.260 

1  FC804 

80 

29.20 

1 . 265 

1 . 304 

1  FC806 

80 

29.20 

1 . 265 

1.304 

FC1007 

100 

29 . 20 

1 . 310 

1.358 

FC1009 

100 

29.20 

1.310 

1.358 
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QUALITATIVE  FLOW  CHARACTERISTICS 


Schlieren  photography  was  used  for  qualitative  flow  analysis. 
Images  of  the  flow  were  captured  on  Polaroid  Type  52  film.  Shock 
patterns  over  the  flat  plate  can  be  seen  from  the  photographs  shown 
in  Figures  17  -  20.  The  boundary  layer  can  be  seen  to  have 
substantial  activity  indicating  that  the  boundary  layer  is 
turbulent  which  is  in  accordance  with  the  heat  transfer  results 
discussed  in  the  later  portion  of  this  chapter.  From  close 
observation  of  the  photographs  some  flow  activity  beyond  the 
boundary  layer  may  be  observed  which  can  be  speculated  to  be  free 
stream  turbulence.  Such  images  are  more  enhanced  in  photographs  of 
flow  when  using  the  turbulence  generator  which  could  be  a  further 
confirmation  that  it  does  increase  free  stream  turbulence.  However, 
such  speculations  can  be  brought  to  further  substance  by  using  high 
speed  photography. 

The  most  interesting  part  of  Schlieren  photography  was  with 
film-cooling  tests.  As  can  be  seen  from  the  comparison  of  Figures 
17  -  20,  the  film  cooled  boundary  layer  may  grow  substantially  in 
thickness  compared  to  the  uncooled  boundary  layer.  As  will  be 
discussed  in  the  film  cooling  section  of  this  chapter,  an  increase 
in  the  film  cooling  flow  beyond  a  certain  point  causes  an  increase 
in  the  thickness  of  the  boundary  layer  and  has  the  adverse  effect 
of  increasing  the  heat  transfer.  High  speed  photography  could  also 
be  useful  in  optimizing  the  method  of  introducing  the  cooling  air. 

During  the  film-cooling  tests  it  was  observed  from  the  heat 
transfer  data  (to  be  discussed  later)  that  an  instantaneous  spike 
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occurred  at  about  1.5  ms  of  the  shock  passage  across  each  heat 
transfer  gage.  Since  the  heat  flux  gages  could  not  be  seen  by  the 
camera  in  the  test  section,  a  new  test  section  with  the  heat  flux 
gages  visible  through  the  windows,  had  to  be  installed  in  an  effort 
to  capture  the  view  of  the  flow  over  the  heat  flux  gages  at  that 
particular  time. 

FLOW  PROPERTIES  BEHIND  THE  SHOCK  WAVE 

The  passage  of  shock  wave  across  a  medium  produces  abrupt 
changes  in  the  properties  of  the  flow  across  the  medium  .  The 
pressure,  temperature  and  density  of  the  flow  behind  the  shock  wave 
increases  whereas  the  velocity  decreases.  The  accuracy  of  the 
measurements  of  these  properties  determine  the  validity  of  the  heat 
transfer  solution. 

The  normal  shock  tube  relationships  given  by  Equations  (3)  to 
(5)  determine  the  theoretical  pressure,  temperature  and  velocity  of 
the  flow  behind  the  shock  wave. 

The  pressure  behind  the  shockwave  was  measured  experimentally 
by  the  front  and  rear  pressure  transducers  located  at  some  distance 
upstream  of  the  flat  plate.  The  temperature  and  velocity 
measurements,  under  the  present  set  up,  were  difficult  to  make 
because  of  the  very  short  duration  of  the  test  runs.  As  such  the 
theoretical  values  of  temperature  and  velocity,  based  on  the 
measured  shock  speed,  v/ere  used  for  analysis  of  all  experimental 
data . 

Flow  properties  like  viscosity  and  thermal  conductivity  were 
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evaluated  at  reference  temperature  calculated  by  using  Equation 
(15).  The  adiabatic  wall  temperatures  were  evaluated  using  the 
theoretical  temperatures  and  Mach  numbers  behind  the  shock  wave, 
based  on  measured  shock  speed.  The  data  presented  as  Appendix  E 
represent  the  flow  properties  evaluated  from  the  above  mentioned 
equations . 

It  may  be  noted  that  a  large  error  between  the  actual  and 
measured  properties  would  pronouncedly  influence  the  heat  transfer 
solution.  A  discussion  to  such  effect  follows  later  in  this 
chapter . 

Real  complications  of  evaluating  the  data  occurred  with  the 
film  cooling  studies.  Since  the  adiabatic  wall  temperature  of  a 
film  cooled  surface  is  different  than  the  adiabatic  wall 
temperature  of  a  surface  with  no  film-cooling.  Equation  16  could 
not  be  used  for  a  film-cooled  case.  The  transient  nature  of  the 
tests  did  not  allow  the  surface  to  reach  an  adiabatic  condition  and 
thus  neither  could  Equation  (16)  be  used  for  theoretical 
calculation  of  the  adiabatic  wall  temperature  (which  would  have 
been  a  misrepresentation)  nor  was  there  any  other  way  of  getting 
around  the  problem.  The  lack  of  knowledge  of  such  an  important 
property  of  the  flow  prevented  the  evaluation  of  properties  like 
film-cooling  effectiveness,  heat  transfer  coefficient,  Nusselt 
number  and  Stanton  number.  The  absence  of  such  an  erudition 
strictly  limited  the  investigation  of  the  data  obtained  during  this 
part  of  the  study.  A  discussion  to  this  effect  follows  in  the 
section  of  heat  transfer  with  film-cooling. 
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TURBULENCE  INTENSITY 


To  acquire  a  true  heat  transfer  solution  to  the  flat  plate 
problem,  knowledge  of  the  free-stream  turbulence  level  and  the 
effect  of  increasing  free-stream  turbulence  on  heat  transfer  was 
required.  This  part  of  the  study  was  a  repetition  of  the  work 
carried  out  by  Rockwell  (1989).  For  comparison  purposes,  the  same 
methodology,  as  adopted  by  Rockwell  (1989)  from  Oldfield  (1984), 
was  employed  to  determine  the  free-stream  turbulence  level  of  the 
flow  in  the  shock  tube. 

Table  V  presents  the  data  obtained  during  the  turbulence  runs 
and  was  used  to  determine  the  turbulence  intensity  of  the  flow. 


TABLE  V:  Turbulence  Data  obtained  from  Turbulence  Runs 


1  Background 

Turbulence 

Turbulence 

Generation  | 

I  Run  No 

Run  No 

Rd(ft) 

^Vris^/'^ieaa^ 

H201 

im 

0.000454 

T101 

6.50 

0.000702 

H202 

0.000455 

T102 

6.75 

0.000508  1 

H203 

0.000545 

T103 

7.00 

0.000593  1 

H204 

7 . 8 

0.000431 

T104 

7 . 25 

0.000468  1 

H205 

7 . 9 

0 . 000403 

T105 

7  50 

0.000524  1 

H206 

8.0 

0 . 000480 

T106 

7.75 

0.000448  1 

J  H207 

8.1 

0.000592 

T107 

8.00 

0.000586  1 

H208 

8 . 2 

0.000572 

T108 

8.25 

0.001250  1 

H209 

8 . 3 

0.000582 

T109 

8.50 

0.001790  1 
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Plots  of  the  ratio  ( )  against  the  overheat  parameter 
are  shown  in  Figures  21  and  22.  Figure  2l  shows  a  quadratic  curve 
fit  of  the  background  turbulence  data  whereas  Figure  22  shows  a 
quadratic  curve  fit  of  the  data  obtained  when  the  turbulence 
generator  was  used.  Table  VI  shows  the  turbulence  intensities 
obtained  from  the  quadratic  curve  fits  of  Figures  21  and  22. 


Table  VI:  Measured  Turbulence  Intensities 


Background  Turbulence 

Turbulence  Injection  | 

9 . 4% 

11.6%  1 

The  background  turbulence  level  of  the  flow  behind  the  shock 
wave  was  found  to  be  a  relatively  high  9.4%.  The  turbulence 
generator  did  not  help  much  in  increasing  the  turbulence  level  of 
the  flow.  The  increased  turbulence  intensity,  due  to  turbulence 
generation,  was  found  to  be  ii.6%.  These  values  of  turbulence 
intensities  correspond  extremely  well  with  the  results  of  Rockwell 
(1989)  and  the  forthcoming  heat  transfer  data. 

It  would  be  worthwhile  to  note  that  though  the  turbulence 
intensity  values  of  this  study  and  the  ones  obtained  by  Rockwell 
(1989)  correlate  extremely  well,  the  verity  of  the  matter  that  the 
hot  film  exhibit  vibrations  after  the  passage  of  a  shock  wave 
thwart  the  confidence  in  the  data.  The  use  of  a  probe  specially 
designed  to  dampen  dynamic  loadings  caused  by  the  passage  of  a 
shock  wave  can  substantially  increase  our  confidence  in  the  data. 
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Such  probes  can  be  provided  by  TSI  corporation  on  special  orders. 


HEAT  TRANSFER  WITHOUT  FILM-COOLING 

The  passage  of  a  shock  wave  across  a  flat  plate  in  the  shock 
tube  initiates  a  high  temperature  flow,  following  the  shock  wave. 
The  resulting  heat  flow  into  the  plate  causes  its  temperature  to 
increase.  This  change  in  the  plate  temperature  is  sensed  by  the 
heat  flux  gages.  The  response  of  the  heat  flux  gages  in  terms  of 
temperature  change  is  illustrated  by  Figures  23  to  25.  It  can  be 
observed  from  these  figures  that  a  step  change  in  the  environment 
produces  a  parabolic  response  by  the  heat  flux  gages.  Though  the 
difference  between  the  temperatures  of  the  flow  and  the  plate  is 
quite  high,  the  rise  in  the  plate  temperature  is  very  nominal 
because  of  the  very  short  duration  of  the  test  runs. 

The  output  of  the  heat  flux  gages,  after  filtering  and 
amplification,  were  converted  into  analog  voltages  by  the  analog 
circuits  and  finally  into  actual  heat  flux  by  using  Equation  (23). 
Figures  26  to  28  represent  the  plots  exhibiting  the  heat  flux  into 
the  flat  plate  as  a  function  of  time.  Note  that  all  these  figures 
show  three  distinct  regions.  The  region  immediately  following  the 
passage  of  the  shock  wave  in  this  impulsively  started  flow 
initially  have  no  boundary  layer  and  therefore  the  rate  of  heat 
transfer  to  the  flat  plate  is  high.  Rapidly  following  the  shock 
wave  is  the  formation  of  an  unsteady  boundary  layer  which  starts 
reducing  the  heat  flux  to  the  plate.  This  unsteady  boundary  layer, 
however,  transitions  fairly  quickly  into  a  steady  boundary  layer 
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which  may  still  be  laminar.  This  whole  region  is  denoted  as  laminar 
region  in  Figure  26.  The  small  bump  in  the  data  within  this  region, 
which  is  consistent  for  all  the  heat  flux  gages,  may  well  be  the 
transitioning  of  the  unsteady  boundary  layer  into  steady  laminar 
boundary  layer.  Following  the  laminar  region  is  the  region  of  low 
heat  transfer  to  the  flat  plate  which  is  denoted  by  transition 
region  in  Figure  26.  The  boundary  layer  within  this  region 
transitions  from  laminar  to  turbulent.  The  third  region,  denoted  by 
turbulent  region,  is  a  region  of  fairly  constant  heat  transfer.  The 
boundary  layer  in  this  region  is  fully  turbulent. 

Figures  23  -  25  are  also  indicative  of  the  advantage  of  using 
the  heat  flux  analog  circuit  in  which  case,  unlike  other  analytical 
techniques,  instantaneous  values  of  heat  flux  can  be  calculated 
without  getting  negative  values  of  heat  flux. 

To  determine  properties  like  Nusselt  number,  and  Stanton 
number  of  the  flow,  an  average  heat  flux  value  for  the  turbulent 
region  was  determined  from  the  heat  flux  versus  time  plots.  The 
Nusselt  number  and  Stanton  number  for  each  run  were  then  determined 
by  using  the  following  equations: 


A.=- 


or. 


(T  -T  ) 

'  av  nr' 


(oo) 


Mi, 


St=-  * 


[i  /  / 


(BJ) 
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From  the  results  of  the  free-stream  turbulence  it  was 


established  that  the  turbulence  intensity  of  the  flow  with  and 
without  turbulence  generation  was  11.6%  and  9.4%  respectively. 

The  presence  of  turbulence  in  the  free-stream  flow  is  expected  to 
affect  the  rate  of  heat  transfer.  Such  an  effect  can  be  seen  in 
Figures  29  and  30. 

Figure  29  shows  Stanton  number  as  a  function  of  Reynolds 
number.  The  solid  and  the  dotted  lines  represent  the  theoretical 
Stanton  number  as  a  function  of  Reynolds  number  with  zero  and  9.4% 
free  stream  turbulence  respectively.  The  solid  line  was  obtained 
from  Equation  46  whereas  the  dotted  line  was  obtained  from  Equation 
33.  The  experimentally  calculated  Stanton  numbers,  represented  as 
data  points  on  the  Figure  29,  correspond  fairly  well  to  the  9.4% 
free-stream  turbulence  curve. 

Figure  30  is  a  representation  of  similar  data  obtained  from 
flow  with  11.6%  free-stream  turbulence.  The  heat  transfer  data,  in 
terms  of  Stanton  number,  seem  to  have  a  fair  agreement  with  the 
theoretical  data. 

Such  a  compatibility  between  experimental  and  theoretical 
results  validates  the  authenticity  of  the  experimental  data.  A 
comparison  of  this  data  with  Rockwell  (1989)  reveals  the  importance 
of  calibration  of  the  heat  flux  gages  for  V  ( oCpk)  .  The  profound 
scatter  in  Rockwell's  data  is  fairly  attributable  to  the  fact  that 
he  used  a  book  value  of  1388.84  J/s*’^-m^-°K  for  all  the  heat 
transfer  gages  whereas  in  actual  the  V  (pCpk^)  value  of  heat  flux 
gages  is  quite  different  than  the  book  value  (refer  to  calibration 
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curves  in  Appendix  B) .  The  use  of  correct  values  of  the  calibration 
constants  reduces  this  scatter  to  a  great  extent  and  can  be  seen 
from  the  comparison  of  Figure  29  with  the  data  of  Rockwell  (1989) 
shown  in  Figure  31. 

HEAT  TRANSFER  WITH  FILM-COOLING 

To  determine  the  effects  of  film-cooling  on  the  heat  transfer 
solution,  cool  air  was  injected  into  the  turbulent  boundary  layer 
formed  on  the  surface  of  the  flat  plate  due  to  the  passage  of  hot 
air  following  the  shock  wave.  Several  test  runs,  at  three  different 
driver  pressures  of  60,  80  and  100  in  of  Hg,  were  made.  Assuming 
very  slight  variations  in  the  properties  of  the  primary  flow  from 
run  to  run  for  each  pressure,  coolant  air  was  blown  into  the 
boundary  layer  at  pressures  such  that  any  desired  blowing  ratio 
could  be  maintained.  The  pressure  of  the  coolant  air  was  maintained 
through  the  film  cooling  system  described  previously.  The  blowing 
ratio  was  determined  from  Equation  33,  whereas  the  properties  of 
both  the  coolant  and  the  primary  flows  were  determined  from  a 
series  of  equations  stated  in  Chapters  II  and  IV. 

Examining  the  output  of  the  heat  flux  gages  in  terms  of  the 
rise  in  temperature  of  the  flat  plate.  Figure  32  shows  the  outputs 
of  Gage  No  2  for  four  different  test  runs  with  three  different 
blowing  ratios  but  the  same  driver  pressure.  With  an  increase  in 
the  blowing  rate,  the  rate  of  change  of  temperature  decreases. 
However,  the  output  curve  with  the  highest  blowing  ratio  shows  that 
the  temperature  of  the  plate  has  started  increasing  though  the  rate 
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of  change  of  temperature  has  decreased.  This  is  demonstrative  of 
the  fact  that  there  comes  a  point  during  the  process  of  film¬ 
cooling  beyond  which  any  further  increase  in  the  blowing  ratio 
starts  showing  adverse  effects.  It  may  be  mentioned  here  that  the 
slight  shift  in  the  initial  values  of  the  data  can  be  attributed  to 
the  slight  changes  in  the  ambient  conditions  taking  place  during 
the  span  of  time  in  which  the  data  was  taken. 

Analyzing  next  the  effects  of  film-cooling  on  heat  transfer 
into  the  flat  plate;  Figures  33  to  42  show  a  few  examples  of  the 
heat  flux  as  a  function  of  time  for  the  film-cooling  runs.  A 
systematic  observation  of  these  graphical  representations  show  that 
the  heat  flux  decreases  with  increasing  blowing  ratio.  However,  as 
stated  before,  a  stage  does  arise  where  any  further  increase  in  the 
blowing  ratio  does  not  decrease  the  heat  flux  any  more.  Instead, 
the  heat  flux  tends  to  start  increasing.  This  reverse  increment  of 
the  heat  flux  could  not  be  explored  deeply  due  to  hardware 
1 imi tat i ons . 

One  important  change  that  can  be  observed  by  comparison  of  the 
heat  flux  output  of  test  runs  without  film-cooling  to  that  of  the 
output  of  test  runs  with  film-cooling  is  an  impulsive  fluctuation 
in  the  heat  flux  at  about  2.0  msec  after  the  shock  passage  over  an 
individual  heat  flux  gage  in  the  later  case.  The  reason  for  such  a 
phenomenon,  however,  is  not  clear  at  present. 

Ignoring  the  above  mentioned  phenomenon,  the  heat  flux  outputs 
of  the  gages  reveal  another  important  and  understandable  fact  about 
the  film-cooling  effectiveness  in  terms  of  heat  transfer.  A 
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comparison  of  Figures  43  and  44,  both  for  different  gage  locations, 
show  that  the  film-cooling  effectiveness  (not  by  definition  but  by 
its  ability  to  reduce  heat  transfer)  of  the  injected  coolant  air 
decreases  along  the  length  of  the  flat  plate.  An  illustration  to 
this  effect  can  be  seen,  though  grimly,  from  the  limited  data 
presented  in  Figure  45.  The  data  limitation  arises  from  the  fact 
that  only  three  of  the  seven  heat  flux  gages  mounted  along  the 
length  of  the  plate  were  available  for  data  collection. 

The  observations  pertaining  to  the  effects  of  film-cooling  on 
the  rate  of  heat  transfer  to  the  flat  plate  cannot  be  concluded 
with  a  forceful  argument  without  its  comparison  to  theoretical  or 
other  experimental  research. 

Heat  transfer,  in  common,  is  mostly  represented  in  terms  of 
heat  transfer  coefficient  and  Stanton/Nusselt  numbers.  To  determine 
the  Stanton/Nusselt  numbers  the  heat  transfer  coefficient  is  a 
required  element.  For  high  speed  flows,  the  heat  transfer 
coefficient  can  be  determined  from  the  equation: 


which  means  that  knowledge  of  heat  flux  qj^ ,  adiabatic  wall 
temperature  Tj^  and  the  wall  temperature  T^  is  imperative. 

What  is  important  to  note  in  this  study  is  that  the  testing 
facility  -the  shock  tube-  is  transient  in  nature  implying  that  the 
short  duration  of  test  time  does  not  allow  the  test  piece,  in  this 
case  the  flat  plate,  to  reach  an  adiabatic  state  and  thus  the 
adiabatic  wall  temperature  cannot  be  determined. 
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Film-cooling  is 
cooling  effectiveness 
equation : 


most  commonly 
which  for  high 


analyzed  in  terms  of  film- 
speed  flows  is  defined  by  the 


»1  = 


'^2 

T^-T, 


(40) 


As  mentioned  earlier,  lack  of  such  critical  information 
blindfolded  this  analysis  from  expansion  into  some  extremely 
important  aspects  of  analysis  thereby  compelling  it  to  remain 
within  the  limits  of  utilizing  the  only  available  data  of  heat  flux 
from  the  heat  flux  gages. 

As  such  Figures  46  shows  a  plot  of  heat  flux  ratio  "q/qj"  as 
a  function  of  blowing  ratio  "B" .  The  three  curves  represent  the 
effects  of  increasing  blowing  ratio  at  different  driver  pressures, 
in  this  case  60,  80  and  100"  Hg.  The  three  important  aspects  of 
film-cooling,  for  the  range  of  conditions  investigated,  revealed  by 
this  figure  are: 

1.  An  increase  in  the  primary  flow  conditions:  i.e,  pressure, 
temperature  and  Mach  number  of  the  primary  flow,,  the  effectiveness 
of  film-cooling  in  reducing  the  heat  transfer  to  the  flat  plate 
decreases . 

2.  A  maximum  decrease  in  the  rate  of  heat  transfer 
(represented  by  the  minima  of  the  curves)  occurs  at  a  point  where 
the  injected  coolant  flow  gets  choked.  Beyond  this  point,  any 
further  increase  in  the  blowing  ratio  starts  showing  adverse 
effects  in  terms  of  increasing  the  heat  flux  to  the  flat  plate. 

3.  The  point  of  minimum  heat  transfer  shifts  upwards  (higher 
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rate  of  heat  transfer)  and  towards  a  lower  blowing  ratio  by 
increasing  the  primary  flow  conditions. 

4.  The  rate  of  decrease  of  heat  flux  to  the  flat  plate  is 
identical  for  all  the  three  primary  flow  conditions  up  to  a  blowing 
ratio  of  unity.  Beyond  this  point  the  data  starts  segregating  on 
the  basis  of  individual  primary  flow  conditions. 

Table  VI  shows  a  consolidated  form  of  the  above  discussion  in 
terms  of  numerical  values: 


TABLErVI  Film  Cooling  Results 


=  60"  Hg 

_ 1 

=  80"  Hg 

=  100"  Hg  j 

R.NO 

B 

R.NO 

B 

R.NO 

B 

q/qo 

FC602 

0.30 

1.00 

FC801 

0.34 

1.00 

FC1001 

0.37 

1.00 

FC604 

0.88 

FC803 

1.39 

0.88 

FC1002 

1.02 

0.93 

FC606 

D 

m 

FC805 

2.97 

0.68 

FC1003 

1 . 99 

0.85 

FC608 

3.93 

0.58 

FC807 

3.23 

0.60 

FC1005 

2 .81 

0.78 

FC6i0 

5.01 

0.52 

FC808 

3.45 

0.  56 

FC1006 

3.08 

0.68 

FC612 

6.63 

0.61 

FC810 

3.90 

0.58 

FC1007 

3.46 

0.70 

FC614 

7,93 

0.64 

FC811 

4.57 

0,60 

FC1008 

3.83 

0.73 

FC616 

9.22 

0.66 

FC812 

5.02 

0.64 

FC1010 

4.58 

0.82 

Note:-  The  shaded  portion  of  the  Table  VI  indicate  that  the  flow  is 
choked . 
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It  may,  however,  be  noted  here  that  the  qualitative  analysis 
of  the  effects  of  film-cooling  on  heat  transfer  from  the  schlieren 
photography  reveal  a  somewhat  different  viewpoint.  A  systematic 
view  of  the  schlieren  photographs  shown  in  Figures  18  -  20  show 
that  by  increasing  the  blowing  ratio  of  the  film-cooling  flow  to 
about  0.8  the  boundary  layer  on  the  flat  plate  increases  with  no 
observable  indication  of  an  increase  in  turbulence  or  the 
separation  of  the  boundary  layer.  Any  further  increase  in  the 
blowing  ratio  beyond  0.8  separates  the  boundary  layer  from  the 
surface  of  the  plate.  Though  the  boundary  layer  quickly  attaches 
back  to  the  surface,  the  turbulence  withen  the  boundary  layer 
increases  significantly.  This  should  indicate  a  decrease  in  the 
effectiveness  of  the  film-cooling  flow.  However  the  quantitative 
analysis  of  the  data  obtained  in  terms  of  heat  transfer  from  the 
heat  flux  gages,  discussed  in  the  preceding  paragraphs,  contradict 
this  qualitative  analysis.  A  better  understanding  of  such 
contradictions  could  be  possible  by  using  high  speed  photography 
for  qualitative  analysis. 

In  view  of  the  data  available  it  can  only  be  speculated  that 
film  cooling  is  most  effective  in  high  speed  flows  at  the  early 
choked  conditions  of  the  film  cooling  flow  through  discrete  holes. 
However,  a  complete  picture  of  the  scenario  can  only  be  determined 
after  obtaining  the  knowledge  of  properties  like  the  adiabatic  wall 
temperature  of  the  plate  and  the  heat  transfer  coefficient  of  the 
flow  with  film  cooling  or  by  determining  empirical  relationships 
which  would  allow  estimation  of  these  properties. 
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VII  CONCLUSIONS 


1.  Calibration  of  the  properties  of  sensing  elements  that  have  a 
direct  effect  on  the  outcome  of  a  solution  is  extremely  necessary 
and  important.  The  effect  of  correct  calibration  constants  could 
seriously  be  felt  to  have  a  positive  influence  on  the  heat  transfer 
solution. 

2.  The  new  method  of  calibrating  the  heat  flux  gages  works 
extremely  well  and  its  results  are  backed  up  strongly  by  the  heat 
transfer  solution. 

3.  The  heat  transfer  solution  of  the  flat  plate,  with 
consideration  for  free-stream  turbulence  and  calibration  constants, 
compares  well  with  the  theoretical  solution. 

4.  Even  though  the  film-cooling  effectiveness  could  not  be 
determined  from  the  present  experimental  set-up;  it  can  be 
established  from  the  data  that  the  heat  transfer  to  a  surface  can 
be  reduced  by  as  much  as  50%  with  film-cooling. 

5.  An  increase  in  the  primary  flow  conditions  results  in  a 
decrease  in  the  effectiveness  of  the  film-cooling  flow  in  reducing 
the  heat  flux  to  the  flat  plate. 

6.  The  point  of  minimum  heat  transfer  due  to  film-cooling  shift 
towards  a  high  heat  transfer  rate  and  lower  blowing  ratio  by 
increasing  the  primary  flow  conditions. 
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VIII  RECOMMENDATIONS 


The  observations  made  during  the  course  of  this  study  lead  to 
the  following  recommendations: 

1.  Further  studies  on  measuring  the  turbulence  level  of  the  flow 
in  the  shock  tube  may  be  carried  out  using  special  purpose  hot 
film/wire  probes  with  a  capacity  to  dampen  the  effects  of 
vibrations.  This  would  relieve  the  suspicion  of  relying  on 
turbulence  data  which  might  be  having  vibration  effects  within  it. 

2.  Since  the  boundary  layer  region  has  the  main  effect  on  the 
heat  transfer  solution  of  the  flat  plate  it  would  be  worthwhile  to 
measure  the  turbulence  level  of  the  boundary  layer  itself. 

3.  In  view  of  the  last  conclusion  stated  in  chapter  VII  it  is 
deemed  necessary  that  high  speed  photography  be  used  for  the 
qualitative  analysis  of  the  film-cooling  flow. 

4  Compare  the  effects  of  film-cooling  on  heat  transfer  by  using 
both  the  90°  and  inclined  film-cooling  injection  holes. 

5.  Find  a  way  to  either  determine  the  adiabatic  wall  temperature 
of  the  flat  plate  in  a  transient  facility  like  the  shock  tube,  or 
a  correlation  that  would  lead  to  establishing  properties  like  the 
Stanton  number  and  Nusselt  number  of  the  flow  with  film  cooling 
affects . 


53 


BIBLIOGRAPHY 


Ammari,  H.d.,  et  al .  "The  effect  of  Density  ratio  on  the  Heat 

Transfer  Coefficient  from  a  Film  Cooled  Flat  Plate",  ASME  Gas 
Turbine  and  Aero  Engine  Congress  and  Exposition,  Toronto, 
Canada,  June  4-8,  1989. 

Blair,  M.F.  "Influence  of  Free  Stream  Turbulence  on  Turbulent 
Boundary  Layer  Heat  Transfer  and  Mean  Profile  Development, 
Part  II  -  Analysis  and  Results.",  Journal  of  Heat  Transfer. 
November,  1983. 

Blasius,  H.  "Grenzschi chten  in  Flussi gkei ten  mit  kleiner  Reibung.", 
Z.  Math  Physics,  56,  1908.  English  Translation  in  NACA  TM 

1256,  1956. 

Bogdan,  Leonard  B.  "High  Temperature,  Thin  Film  Resistance 
Thermometers  For  Heat  Transfer  Measurements”,  Cornel  1 
Aeronautical  Laboratory  Inc.  June,  1967. 

Bradshaw,  P.  "An  Introduction  to  Turbulence  and  Its  Measurement", 
Perqammon  Press,  1971. 

Chapmann,  Alan  J.  and  W.F.  Walker.  "Introductory  Gas  Dynamics", 
Holt,  Rinehart  and  Wilson  Inc,  New  York.  1971. 

Crawford,  M.E,  et  al,  "Film  Cooling  Effectiveness  Downstream  of  a 
Single  Row  of  Holes  with  Variable  Density  Ratio",  Gas  Turbine 
and  Aeroengines  Congress  and  Exposition,  June  11-14,  1990. 

Davis,  W.R.  and  L.  Bernstein,  "Heat  Transfer  and  Transition  to 
Turbulence  in  a  Shock  Induced  Boundary  Layer  on  a  Semi 
Infinite  Flat  Plate",  Journal  of  Fluid  Mechanics,  Volume  36: 
Part  1.  1969. 

Dillon,  R.e.  and  H.T.  Nagamatsu.  "Heat  Transfer  and  Transition 
Mechanism  on  a  Shock  Tube  Wall",  AIAA  Journal,  22:  no  ll, 
November,  1984. 

Dunn,  M.G.  and  R.J.  Goldstein.  "Measurements  in  Heat  Transfer  to 
Gas  Turbine  Components",  Technical  Report  AFAPL-TR-77-66 
Calspan  Corporation,  New  York.  October,  1977. 

Eckert,  E.P.G  and  R.J.  Goldstein.  "Measurements  in  Heat  Transfer", 
Second  Edition,  McGraw  Hill  Book  Company.  NY.  1976. 

Gaydon,  A.G.  and  I.R.  Hurle.  "Shock  Tube  in  High  Temperature 
Chemical  Physics",  Reinhold  Publishing  Corporation.  NY.  1963. 

Glass,  I. I.  "Shock  Tubes,  Part  I:  Theory  and  Performance  of  Simple 
Shock  Tubes",  UTIA  Review  No  12,  Toronto . Canada :  Institute  of 


54 


Aerophysics,  University  of  Toronto.  1958. 


Goldstein,  R.J.,  et  al .  "Advances  in  Heat  Transfer",  Academic  Press 
NYr  1971. 

Goldstein,  R.J.,  et  al .  "Film  Cooling  Following  Injection  Through 
Inclined  Circular  Holes",  NASA  Report  No  CR-72612.  November, 
1969 . 

Han,  J.C.  and  A.B.  Mehendale.  "Influence  of  High  Mainstream 
Turbulence  on  Leading  Edge  Film  Cooling  Heat  Transfer",  Gas 
Turbine  and  Aeroengines  Congress  and  Exposition.  June  11-14, 
1990. 

Hartunian,R. A. ,  et  al  "Boundary  Layer  Transition  and  Heat  Transfer 

in  Shock  Tubes" . Journal  of  Aerospace  Sci ences , 27 : 587 .  1960 

Hinze,  J.O.  "Turbulence:  An  Introduction  to  Its  Mechanism",  McGraw- 
Hill  Book  Company.  1980. 

Jurgelewicz,  Scott  a.  "Investigation  of  Heat  Transfer  with  Film 
Cooling  to  a  flat  plate  in  a  Shock  Tube",  Masters  Thesis. 
School  of  Engineering,  Air  Force  Institute  of  Technology  (AU) , 
Wright  Patterson  Air  Force  Base,  OH.  December,  1989. 

Kays,  W.M.  and  M.E.  Crawford.  "Convection  Heat  and  Mass  Transfer", 
McGraw-Hill  Book  Company,  NY,  1980. 

Kakac,  s.  et  al .  "Hand  Book  of  Single  Phase  Convective  Heat 
Transfer",  John  Willy  and  sons,  NY.  1987. 

Ko.Y.S.  and  D.Y.  Liu.  "Experimental  Investigation  on  The 
Effectiveness,  Heat  Transfer  Coefficient  and  Turbulence  of 
Film  Cooling",  AIAA.  August,  1980. 

Lam,  H.  and  L.Crocco.  "Shock  Induced  Unsteady  Laminar  Compressible 

Boundary  Layer  on  a  Semi  Infinite  Flat  Plate,  "Princeton  University 

Report  No  428,  1958. 

MacMullin,  R.,  et  al  .  "Free  Stream  Turbulence  From  a  Circular  Wall 
Jet  on  a  Flat  Plate  Heat  Transfer  and  Boundary  Layer",  ASME . 
January,  1989. 

McQueen,  Stephen  M.  "Velocity  and  Transient  Measurements  in  a  Shock 
Tube  Using  a  Hot  Wire  Anemometer",  Masters  Thesis.  School  of 
Engineering,  Force  Institute  of  Technology  (AU),  Wright 
Patterson  Air  Force  Base  OH,  December,  1984. 

Mirels,  J.  "Boundary  Layer  Behind  Shock  of  Thin  Expansion  Wave 
Moving  into  a  Stationary  Fluid",  NACA  TN  3712,  1956. 

Novak,  Joseph  T.  "Investigation  f  Heat  Transfer  To  a  Flat  Plate  in 


55 


a  Shock  Tube",  Masters  Thesis,  School  of  Engineering,  Air 
Force  Institute  of  Technology  (AU),  Wright  Patterson  Air  Force 
Base,  OH.  December,  1987. 

Oldfield,  M.L.G.,  et  al  .  "On  Line  Computer  for  Transient  Turbine 
Cascade  Instrumentation",  IEEE  Transactions  on  Aerospace  and 
Electronic  Systems  Vol  AES-14.  No  5.  September,  1978. 

Prandtl.L.  "Uber  Flussigkeitsbewegung  Bei  Sher  Kleiner  Reibung 

Proc",  Third  International  Math  Congress.  Heidelberg:  484-491,  1904 

Rockwell,  Richard  K.  "Transient  Heat  Transfer  Measurements  on  a 
Flat  PLate  in  Turbulent  Flow  Using  an  Electrical  Analog", 
Masters  Thesis.  School  of  Engineering,  Air  Force  Institute  of 
Technology  (AU) ,  Wright  Patterson  Air  Force  Base,  OH. 
December,  1989. 

Schlichting,  Hermann.  "Boundary  Layer  Theory  (Seventh  Edition)", 
McGraw-Hill  Book  Company.  1989. 

Schmitz,  L.S.  "Nonlinear  Analog  Network  to  Convert  Surface 
Temperature  to  Heat  Transfer",  Cornell  Aeronautical 
Laboratory.  Inc,  Buffalo.  NY.  June,  1963. 

Schultz,  D.L.  and  T.V.  Jones.  "Heat  transfer  Measurements  in  Short 
Duration  Hypersonic  Facilities",  AGARDograph  No.  165. 
February,  1973. 

Shapiro,  Ascher  H.  "The  Dynamics  and  Thermodynamics  of  Compressible 
Fluid  Flow,  Vol.  2,  The  Ronald  press  company.  NY,  1987. 

Simonich,  J.C.  and  P.  Bradshaw.  "Effect  of  Free  Stream  Turbulence 
on  Heat  Transfer  Through  Turbulent  Boundary  Layer",  Journal  of 
Heat  Transfer.  November,  1978. 

Skinner,  George  t.  "Analog  Network  To  Convert  Surface  Temperature 
to  Heat  Flux",  Report  no  CAL-100.  Cornell  Aeronautical 
Laboratory , Inc ,  Buffalo,  NY.  February,  1960. 

Smith,  Bret  J.  "Investigation  of  Heat  Transfer  to  a  Sharp  Leading 
Edged  Flat  Plate  Using  a  Shock  Tube",  Masters  Thesis,  School 
of  Engineering,  Air  Force  Institute  of  Technology.  Wright 
Patterson  Air  Force  Base  OH.  December,  1986. 

Spence,  D.A.  "Velocity  and  Enthalpy  Distribution  in  Compressible 

Turbulent  Boundary  Layer  on  a  Flat  Plate",  Journal  of  Fluid 

Mechanics ,18:117,  1960. 

Zucrow,  Maurice  J.  and  J.D.  Hoffman.  "Gas  Dynamics,  Vol.  I",  John 
Willy  and  Sons ,  NY,  1976. 


56 


TEMPERATURE  T 


TIME 


(a)  Flow  pattern  in  a  shock  tube  at  t=t| 

(b)  Pressure  distribution  in  a  shock  tube  at  t=t2 

(c)  Temperature  distribution  in  a  shock  tube  at  t=t| 

(d)  Physical  plane,  showing  shock  wave,  rarefaction  wave,  and  contact  discontinuity 
(Shapiro,  19S7:  1007) 

Figure  2:  Illustration  of  Flow  Properties/Pattems  in  a  Shock  Tube 
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;  noiindary  Layer  Formalion  Hcliind  a  Shock  Wave 
(SchllclUing,  1987) 
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Fiy.  6:  Film  C(M)liny  and  riirhulcucc  (Jcneratlon  Syslc'ins 


MIRROR  1 


Fig.  7:  Sclilieren  System  Confiyiiralion 


Fig,  8;  Pressure  Transducers  and  Heat  Flux  (Jage  Locations  in  the  Shock  Tube 
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Fig.  9:  Schematic  of  Heat  Flux  Analog  Circuit 


Fig.  10:  Bridge  Circuit  for  Calibration  and  Collection  of  Data  from  H-F  Gages 
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Fig.  12:  Calibration  Circuit  to  Determine  ^pC^k  of  Heat  Flux  Gages 
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Fig,  13:  Schematic  of  Calibration  Circuit  to  Determine  ^pCp/t  of  Heat  Flux  Gages 
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Figu'^e  14:  Parabolic  Output  of  Heat  Flux  Gage  No  4 
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Fig'jre  15:  Output  o'  TS!  1214-11.5  Hotwre 


Figure  17:  schlieren  Photograph  of  Flow  Behind  a  Shock  Wave 
:No  Film  Cooling  or  Turbulence  Generation 


Figure  18:  Schlieren  Photograph  of  Flow  Behind  a  Shock  Wave 

;With  Film  Cooling  B  =  0.20 
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Figure  l9;Schlieren  Photograph  of  Flow  Behind  a  Shock  Wave 

:With  Film  Cooling  B  =  0.80 


Figure  20:  Schlieren  Photograph  of  Flow  Behind  a  Shock  Wave 

:With  Film  Cooling  B  =  2.00 
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Figure  21:  Hot  Film  Powc'  Output  as  a  Function  o"?^ 
Ovcheat  Pc'ameter  (NO  Turbulence  injection) 
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Figure  24:  Temperature  change  as  a  function  of  time 

Gage  No  4 
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Figure  25:  Temperature  change  as  a  function  of  time 

Gage  No  6 
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Figure  26:  Heat  F!ux  as  c  Function  of  Time 
Run  No  R002 
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Figure  27:  Heat  Flux  as  a  Function  of  Time 
Run  No  R005 
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28:  Heat  Fiux  as  a  Function 
Run  No  R007 
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Figure  29:  Stanton  Numbe*-  as  a  Eunct’on  of  Reynolds  Mumbc'- 

No  Turbulence  Generation 
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Igure  30:  Stanton  Number  as  c  ^unction  of  Reynolds  Numbe 
With  Turbulence  Generation 
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Figure  31:  Stanton  Number  os  o  Function  Reynolds 
Number,  No  Turbulence  Generation 

(Rockwell,  1989) 
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Figure  32:  Temperature  Change  as  a  Function  of  Time 

with  Mim  Cooling 


85 


GAGE  NO  2  : 

Turbulence  intensity  9.6s5  : 

Blowing  Ratio  0.2  : 

I  I  I  1  ]  I  1  ]  1  i  i  1  i  I  I  -Li  i  .1  1  J  I  J  1  I  I- 

6.0  8.0  10. 
TIME  (msec) 


Figure  33;  Heat  Flux  as  a  Function  of  Time:  FC801 
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Figure  34:  Heat  Flux  as  a  Function  of  Time;  FC802 
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Figure  To:  Heat  Flux  as  a  Function  of  Time:  FC803 
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Flux  as  a  Function  of  Time:  FC804 
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Figure  38:  Heat  Flux  as  a  Function  of  Time:  FC806 
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Figure  39:  Heat  Flux  as  a  Function  of  Time:  FC807 


Figure  40:  Heat  Flux  as  a  Function  of  Time:  FC808 
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Figure  41:  Heat  Flux  as  a  Function  of  Time:  FC811 
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Figure  43:  Heat  Flux  as  a  Function  of  Time:  Gage  No  2 
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Figure  46:  Heat  Flux  Ratio  as  a  Function  of  Blowing  Ratic 


APPENDIX  A 


Derivation  for  i/oCpkpj^j.^ 

One  dimensional  heat  conduction  equation  is  expressed  as: 

(A  1) 

at  pc^dx‘ 

where  T  is  the  difference  between  instantaneous  and  initial 
temperature,  i.e  T(x,0)  =  0. 

For  the  case  of  a  step  function  in  heat  transfer  rate  the 
boundry  conditions  are: 

T(x,0)  =  q(x,0)  =0 

(A.  2) 

q(0,t)  =  q^  =  constant 

The  solution  to  equation  A-1  can  be  found  in  many  references 
as : 


r(o,  t)  = 


2cj 

y/n(pCj^) 


(A. 3) 


Let  Rj,  be  the  initial  resistance  of  the  platinum  foil  in  the 
heat  flux  gage.  By  passing  a  current  of  known  amplitude,  the 
resistance  of  the  platinum  foil  changes  to  R|. 

R{  and  Rg  are  related  to  each  other  by  the  relationship: 

R{  =  Rg  (1  +  aT)  (A. 4) 

which  would  imply  that: 


T= 


AR 

aR^ 


(A. 5) 


aR  can  be  related  to  the  voltage  signal  across  the  heat  flux 
gage  AV(t),  and  the  constant  current  through  the  film  i,  by: 

A. 


aR  -  aV(  t )  /i  (A,  6) 

Substituting  Equations  A. 4  and  A. 5  into  Equation  A. 2  one  obtains 


after  some  manipulation  the  following  equation: 


A  o- 


(A. 7) 


As  explained  in  the  calibration  procedure  for  ^  pCpk ,  the  heat 
flux  gage  is  first  calibrated  in  air  and  the  factor  aRj/ZT  is 
obtained  from  the  resultant  curve.  If  the  gage  is  then  calibrated 
with  a  small  droplet  of  glycerine  on  top  of  the  heat  flux  gage,  the 
heat  produced  due  to  the  passage  of  current  would  flow  both  into 
the  glycerine  and  the  gage  substrate.  If  a  fraction  "m"  of  the  heat 
generated  diffuses  into  glycerine  and  is  the  slope  of  the 
resultant  curve  then: 


v/T 


(A. 8) 


But  also: 


^PCpk 


subs  _ 


2  (l-/n) 


R^a 


>rt 


(A. 9) 


subs  _  2X  R^ 


R^a 


(A. 10) 


From  equations  A-8,  A-9,  and  A-10,  by  manipulation  we  can 


A. 2 


obtain  the  correct  relationship  for  the  bulk  thermal  diffusivity  of 
the  substrate  as: 


A. 3 


APPENDIX  B 


CALIBRATION  OF  HEAT  FLUX  GAGES  FOR  VoCpK 

The  calibration  procedure  adopted  to  find  the  bulk  thermal 
diffusivity  of  the  heat  flux  gages  is  as  follows: 

A  bridge  circuitry  shown  in  Figure  12  was  constructed  using 
precision  resistors.  The  heat  flux  gage  under  calibration  was  used 
as  one  of  the  legs  of  the  bridge  circuitry  whereas  a  variable 
resistor  was  used  parallel  to  the  leg  adjacent  to  the  heat  transfer 
gage  so  as  to  ensure  proper  balancing  of  the  bridge. 

A  500n  resistor  was  used  in  series  with  the  bridge  to  keep  the 
current  at  a  nominal  value  of  less  than  2.5  mamp  to  avoid  any 
excessive  passage  of  current  through  the  circuit  resulting  in 
burning  out  the  heat  flux  gage. 

A  5  volt  DC  power  supply  from  a  Tektronix  work  station  was 
used  to  balance  the  bridge.  But  to  pass  a  single  pulse  of  known 
amplitude  through  the  circuit  a  Wavetek  model  signal  generator  was 
used . 

The  output  of  the  bridge  circuitry  was  connected  to  a 
differential  amplifier  for  filtering  and  amplification.  Having 
filtered  the  output  at  10  Khz  and  amplified  with  a  gain  of  250,  it 
was  fed  into  a  digital  oscilloscope  for  display  and  DL-1200  for 
recording.  A  schematic  of  the  instrument  interfacing  is  shown  in 
Figure  13. 

The  first  part  of  the  calibration  was  accomplished  by  passing 
a  single  square  pulse  of  5  volts  with  a  period  and  width  of  10  msec 
and  5  msec  respectively  through  the  bridge  circuitry.  This  allowed 

to  raise  the  temperature  of  the  heat  flux  gage  resulting  in  a 

B .  : 


change  in  resistance  of  the  platinum  foil  of  the  gage  which  in  turn 
unbalanced  the  previously  balanced  bridge. 

The  amplified  output  of  the  bridge  was  displayed  and  recorded 
simultaneously  on  the  digital  oscilloscope  and  the  DL-1200 
respectively.  Figure  B.l  shows  the  output  voltage  of  the  heat  flux 
gage  w.r.t  time. 

The  second  step  of  the  calibration  involved  a  repetition  of 
the  same  procedure  except  that  the  heat  flux  gage  was  covered  with 
a  droplet  of  glycerine  whose  bulk  thermal  diffusivity  was  known. 

Passage  of  the  current,  this  time,  not  only  allowed  the 
generated  heat  to  flow  into  the  substrate  of  the  gage  but  also  into 
glycerine.  The  output  of  the  heat  flux  gage  was  once  again  recorded 
and  displayed  on  the  DL-1200  and  the  digital  oscilloscope 
respectively.  Figure  B.2  shows  the  o..tput  voltage  of  the  heat  flux 
gage,  covered  with  glycerii.a,  w.r.t  time. 

Figures  B.l  and  B..:  shov;  tha^-  the  output  of  the  heat  flux  gage 
in  both  the  experiments  is  parabolic  with  time.  This  would  imply 
that  a  plot  of  the  output  voltage  w.r.t  fT  would  be  a  linear  curve. 
Figure  B.3  illustrate  an  example  of  such  a  plot. 

The  slope  of  the  curves  illustrated  in  figure  B.3,  i.e  aVj//T 
for  the  experiment  with  heat  flux  gage  in  air  environment  and 
aVj/iTT  for  the  gage  in  glycerine  environment,  are  used  in  Equation 
36  to  find  out  the  correct  value  of  the  bulk  thermal  diffusivity  of 
the  heat  flux  gages.  The  calibration  curves  for  each  gage  are  shown 
in  Figures  B.3  to  B.5. 
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Figure  B.1  :  Parabolic  Output  of  Heat  Flux  Gage  No  4 
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Figure  B.2:  Parabolic  Output  of  Heat  F!ux  Gage  No  4 
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Figure  B.5:  Calibration  Curves  for  Gage  No  6 
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APPENDIX  C 


CALIBRATION  OF  HEAT  FLUX  GAGES  FOR  TEMPERATURE  COEFFICIENT 

The  heat  flux  gages  were  calibrated  for  their  temperature 
coefficient  following  the  procedure  mentioned  in  section  IV. 

The  output  voltage  of  the  bridge  circuitry  was  plotted  against 
temperature.  The  calibration  curves  are  shown  as  Figures  C.l  to 
C.3. 

The  linear  relationships  obtained  as  a  result  of  fitting  a 
curve  through  the  data  are: 

HEAT  FLUX  GAGE  NO  2 

V  =  0.0241806*T  -  1.501 

HEAT  FLUX  GAGE  NO  4 

V  =  0.0236467*T  -  1.6601 

HEAT  FLUX  GAGE  NO  6 


V 


0.01788V*T 


1.42726 
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APPENDIX  D 


CALIBRATION  OF  PRESSURE  TRANSDUCERS 
The  pressure  transducers  used  during  this  study  were 
calibrated  using  an  AMETEK  Model  HK-500  Pneumatic  Dead  Weight 
Tester.  The  output  voltage  of  the  pressure  transducers  was  plotted 
as  a  function  of  pressure  (psig). 

The  calibration  curves  are  plotted  un  Figures  D.l  to  D.3.  The 
least  square  fit  of  the  calibration  curves  resulted  in  the 
following  equations: 

Forward  Pressure  Transducer 
P  =  19.887V  -  0.25537 

Rear  Pressure  Transducer 
P  =  29.504V  +  0.20712 

Film  Cooling  Flow  Pressure  Transducer 
P  =  32.465V  -  0.05925 
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Figure  D.2:  Calibration  Curve  For  Rear  Pressure  Transducer 
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Figure  D.3:  Calibration  Curve  For  F.C  Pressure  Transducer 
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APPENDIX  E 


NUMERICAL  SUMMARY  OF  HEAT  TRANSFER  RUNS 


Us  i  ng 
of  some  of 


the  applicable  equations  stated  in  the  text,  a  summary 
the  Heat  Transfer  tests  is  given  in  the  next  few  pages. 


E.l 


RUN  NO  R002 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 

Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


60.00 

inches  of  Hg  (Gauge) 

29.07 

inches  of  Hg 

44.69 

inches  of  Hg 

23.00 

degrees  C 

62.48 

degrees  C 

344.87 

meters/sec 

416.80 

1.209 

1.286 

meters/sec 

109.54 

0.298 

meters/sec 

67.74 

degrees  C 

44.43 

degrees  C 

1.634 

Kg/m^ 

19.34E- 

6  Pa .  s 

1005.5 

J/Kg  °K 

3  W/m’  °K 

27.45E- 

0.8824 

0.7087 

RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m^)  1 

0.06112 

0.5656 

0.00233 

933.68 

419.35 

18.76  1 

0.07104 

0.6575 

0.00226 

1053.11 

406.92 

18.21 

0.09128 

0.8448 

0.00215 

1287.00 

387.02 

17.32 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


X(m) 

Re(  10®) 

St 

Nu 

h 

q(Kw/m^)  1 

0.06112 

0.5654 

0.00335 

1343.70 

603 . 49 

27.00  1 

0.07104 

0.6572 

0.00323 

1504.00 

581 . 14 

26.00 

0.09128 

0.84442 

0.00320 

1915.00 

575.90 

25.77 

E.  2 


RUN  NO  R003 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 

Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


70.00 

inches  of  Hg  (Gauge) 

29.07 

inches  of  Hg 

46.31 

inches  of  Hg 

23.00 

degrees  C 

66.12 

degrees  C 

344.87 

meters/sec 

423.56 

meters/sec 

1.228 

1 . 297 

118.97 

meters/sec 

0.322 

72.33 

degrees  C 

45.80 

degrees  C 

1 . 675 

Kg/m^ 

19 . 38E-6 

Pa .  s 

1005.6 

J/Ka  «K 

27 . 5E-3 

W/m^  °K 

0.8827 

0.7084 

RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m2)  1 

0.06112 

0.6285 

0.00228 

1015.47 

457.03 

22.55  I 

1  0.07104 

0.7305 

0.00221 

1145.36 

443.37 

21.87  1 

1  0.09128 

0.9386 

0.00211 

1399.70 

421 . 69 

20.80  1 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


1  X(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m^)  1 

0.06112 

0.6287 

0.00354 

1576.00 

709.51 

35.00  1 

0.07104 

0.7308 

0.00364 

1937.60 

729 . 78 

36.00  1 

0.09128 

0.9390 

0.00325 

2161 . 77 

651 . 28 

32.13  1 

E.  3 


RUN  NO  R004 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 


=  P4  =  80 
=  Pj  =  29.09 
=  P2  =  49.78 
=  Tj  =  23.00 
=  Tj  =  73.73 
=  aj  =  344.87 
=  Uj  =  437.66 
=  Mj  =  1.269 
=  M^  =  1.323 


inches  of  Hg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  C 
degrees  C 
meters/sec 
meters/sec 


Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


U2  =  138.26 
M2  =  0.370 

Tj^=  82.11 

T  =  50.04 


meters/sec 

degrees  C 
degrees  C 


Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =  1.701  Kg/m^ 

=  U  =  19.57E-6  Pa.S 
=  C  =  1006.00  J/Kq  °K 

=  k  =  27.8E-3  W/m^  °K 

=  r  =  0.8829 

=  Pr  =  0.7081 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m^)  1 

0.06112 

0.7604 

0.00220 

1182.00 

537.83 

31.80 

0.07104 

0.8839 

0.00213 

1333.65 

521.88 

30.85 

0.09128 

1  .  1356 

0.00203 

1629.8 

496 . 37 

29.34 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

Nu 

h 

q(Kw/m^) 

0.06112 

0.7605 

0.00311 

1673.00 

761 . 29 

45.00 

0.07104 

0.8841 

0.00318 

1988.70 

778.21 

46.00 

0.09128 

1  .  1358 

0.00320 

2574.00 

783.80 

46.33 

E.  4 


RUN  NO  R005 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


=  Pj  =  90.00 
=  Pj  =  29.07 
=  Pj  =  52.60 
=  Tj  =  23.00 
=  T2  =  79.78 
=  aj  =  344.87 
=  Uj  =  448.82 
=  Mj  =  1.302 
=  Mj  =  1.347 
=  Uj  =  153.19 
=  M,  =  0.407 
=  T,„=  90.10 
=  53.45 


inches  of  Hg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  C 
degrees  C 
meters/sec 
meters/sec 


meters/sec 

degrees  C 
degrees  C 


Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =  1.829  Kg/m^ 

=  u  =  19.73E-6  Pa.s 
=  Cp  =  1006.34  J/Kg  °K 

=  k^  =  28.02E-3  w/m"  °K 

=  r  =  0.8831 

=  Pr  =  0.7077 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m^) 

0.06112 

0.8679 

0.00214 

1314.00 

602.83 

40.45 

0.07104 

1.0089 

0.00208 

1482.00 

584.96 

39 . 25 

0.09128 

1 . 2928 

0.00197 

1811.00 

556 . 36 

37.33 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

Nu 

h 

q(Kw/m^) 

0.06112 

0.8684 

0.00307 

1884.11 

864 . 38 

58.00 

0.07104 

1.0094 

0.00296 

2114.45 

834.58 

56.00 

0.09128 

1 . 2970 

0.00292 

2680.22 

823 . 33 

55.25  1 

E.  5 


RUN  NO  R006 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


=  P|  =  100.00 
=  Pj  =  29.07 
=  P2  =  53.61 
=  Tj  =  23.00 
=  T,  =  81.91 
=  aj  =  344.87 
=  Uj  =  455.53 
=  Mj  =  1.32 
=  M^  =  1.37 
=  U2  =  161.63 
=  M2  =  0.4  28 
=  T.„=  93.39 
=  tJ  =  54.75 


inches  of  Hg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  c 
degrees  C 
meters/sec 
meters/sec 


meters/sec 

degrees  C 
degrees  C 


Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =  1.868  Kg/m^ 

=  u  =  19.80E-6  Pa.S 
=  Cp  =  1006.5  J/Kg  °K 

=  k^  =  28. 13E-3  w/m"  °K 

=  r  =  0.8831 

=  Pr  =  0.7074 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10*) 

St 

NU 

h 

q(Kw/m2) 

0.06112 

0.9320 

0.00211 

1390.60 

640.04 

45.05 

0.07104 

1.0833 

0.00205 

1568.47 

621.06 

43,72 

0.09128 

1 . 3919 

0.00195 

1916.76 

590.69 

41.58 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


X{m) 

Re(  10*) 

St 

NU 

h 

qiKw/m^) 

0.06112 

0.9333 

0 . 00304 

1913.70 

923 . 43 

65.00 

0.07104 

1.0841 

0.00280 

2147.32 

850.27 

59.00 

0.09128 

1.3939 

0.00298 

2938 . 56 

905 . 58 

62.80 

RUN  NO  R007 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 

Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  P4  = 
=  P,  = 
=  P2  = 

=  Tj  = 
=  Tj  = 

=  ai  = 
=  = 

=  M,  = 

=  Mj  = 
=  = 

: 

=  T,  = 

=  p  = 

=  u  = 
= 

=  r  = 
=  Pr  = 


110. .00 

inches  ofHg(Gauge) 

29.07 

inches  of  Hg 

52.43 

inches  of  Hg 

23.00 

degrees  C 

87.78 

degrees  C 

344.87 

meters/sec 

463.53 

meters/sec 

1.344 

1 . 390 

172.46 

meters/sec 

0.453 

100.85 

degrees  C 

58.00 

degrees  C 

1.918 

Kg/m^ 

19. 93E- 

6  Pa .  s 

1006.8 

J/Kg  °K 

28. 37E- 

3  W/m'  °K 

0.8833 

0.7072 

RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re( 10S 

St 

NU 

h 

q(Kw/m2)  j 

0.06112 

1 

.0144 

0.00207 

1487 . 85 

690.65 

53.767  1 

0  tci T  ■'  04 

1 

.  1791 

0.00201 

1678.17 

670.175 

52.17  1 

1  0.09128 

1 

.5150 

0.00191 

2050.82 

637 . 405 

49.62 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


X(m) 

Re(  10^) 

St 

Nu 

h 

q(Kw/m^)  1 

0.06112 

1 

.0145 

0.00301 

1605.00 

1001.93 

78.00 

0.07104 

1 

.1792 

0.00270 

2251.60 

899 . 17 

70.00 

0.09128 

1 

.5151 

0.00275 

2946 . 60 

915.82 

71.30 

RUN  NO  R008 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 

Dens i ty 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


120.00 

inches  of  Hg( Gauge) 

29.07 

inches  of  Hg 

57 . 86 

inches  of  Hg 

23.00 

degrees  C 

90.74 

degrees  C 

344.87 

meters/sec 

468.94 

meters/sec 

1 . 360 

1.410 

179.43 

meters/sec 

0.469 

104.88 

degrees  C 

59.70 

degrees  C 

1.951 

Kg/m^ 

20.02E- 

6  Pa .  s 

1007.00 

J/Kg  °K 

28.51E- 

3  W/m^  °K 

0 . 8834 

0.7070 

RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


1  x(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m^)  1 

0.06112 

1 .0721 

0.00205 

1555.00 

725.35 

59.39  1 

0.07104 

1 . 2462 

0.00199 

1754.00 

703.85 

57.63  1 

1  0.09128 

1 . 6012 

0.00189 

2143.00 

669 .43 

54.81  1 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


X(m) 

Re(  10^) 

St 

Nu 

h 

q(KVi/mh  1 

0.06112 

1.0689 

0.00274 

2068.30 

964 . 83 

79.00  I 

0.07104 

1 . 2425 

0.00291 

2556.30 

1025 . 89 

84.00  1 

0.09128 

1 . 5964 

0.00285 

3216 . 70 

1004 . 70 

82.30  1 

E.  9 


RUN  NO  R102 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


=  P^  =  60.00 
=  Pj  =  29.21 
=  P2  =  44.48 
=  Tj  =  24.00 
=  T2  =  63.12 
=  aj  =  345.47 
=  Uj  =  416.34 
=  Mj  =  1.205 
=  Mj  =  1.270 

=  U2  =  108.10 

=  M2  =  0.294 
=  T3^=  68.10 
=  Tj  =  45.60 


inches  of  Hg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  c 
degrees  C 
meters/sec 
meters/sec 


meters/sec 

degrees  C 
degrees  C 


Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =  1.623  Kg/m^ 

=  U  =  19.34E-6  Pa.s 
=  CL  =  1005.5  J/Kg  °K 

=  k  =  27.46E-3  W/m'  °K 

=  r  =  0.883 

=  Pr  =  0.709 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re( 10^ 

St 

NU 

h 

q(Kw/m^)  1 

0.06112 

0.5545 

0.00234 

919.20 

413.00 

18.21  1 

0.07104 

0.6446 

0.00227 

1037.00 

400.84 

17 . 68 

0.09128 

0.8282 

0.00216 

1267.00 

381 . 16 

16 .81 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


X(m) 

Re( 10^ 

St 

NU 

h 

q(KW/m^)  1 

0.06112 

0.5538 

0.00380 

1492 

670.5 

29.6 

0.07104 

0.6439 

0.00342 

1552 

600.0 

26 . 5 

0.09128 

0.8273 

0.00345 

2023 

608 . 73 

26.9  j 

E.  10 


RUN  NO  R103 


Driver  Pressure  =  = 
Driven  Pressure  =  Pj  = 
Pressure  Behind  the  Shock  =  P2  = 
Driver/Driven  Temperature  =  Tj  = 
Temperature  behind  shock  =  T2  = 
Sonic  speed  =  aj  = 
Measured  shock  speed  =  Uj  = 
Measured  shock  Mach  No  =  Mj  = 
Theoretical  shock  Mach  No  =  M^  = 
Flow  velocity  behind  shock  =  U2  = 
Flow  Mach  No  behind  shock  =  M2  = 
Adiabatic  wall  Temperature  =  Tj^= 
Reference  Temperature  =  Tj  = 

Density  =  p  = 
Dynamic  Viscosity  =  p  = 
Specific  Heat  =  Cp  = 
Thermal  conductivity  =  k^  = 
Recovery  Factor  =  r  = 
Prandtl  Number  =  Pr  = 


70.00 

inches  of  Hg  (Gauge 

29 .21 

inches  of  Hg 

46 . 60 

inches  of  Hg 

24.00 

degrees  C 

66 .51 

degrees  C 

345.47 

meters/sec 

424.90 

meters/sec 

1 . 232 

1 . 305 

119.45 

meters/sec 

0.324 

73.60 

degrees  C 

46.80 

degrees  C 

1 . 69 

Kg/m^ 

19. 39E-6 

Pa .  s 

1005 . 6 

27 . 5E-3 

0  0 

0.883 

0.708 

RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m^) 

0.06112 

0.6363 

0.0023 

1025 . 30 

461 . 33 

22.9 

0.07104 

0.7396 

0.0022 

1156 . 40 

447 . 64 

22 . 2 

0.09128 

0.9503 

0.0021 

1413. 18 

425 .75 

21.2  1 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

Nu 

h 

q(Kw/m^)  1 

0.06112 

0.6372 

0.0035 

1578.9 

710.43 

35.2  1 

0.07104 

0.7407 

0.0036 

1887.9 

730.80 

36.3  1 

j  0.09128 

0.9517 

0.0034 

2291.0 

690.20 

34.2  1 

E.  11 


RUN  NO  R104 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


=  P^  =80.00 
=  Pj  =  29.21 
=  P2  =  50.22 
=  Tj  =  24.00 
=  T2  =  74.38 
=  aj  =  345.47 
=  Uj  =  438.20 
=  M.  =  1.272 
=  =  1.330 

=  U2  =  139.2  2 
=  M2  =  0.373 
=  T,„=  83.48 
=  Tj  =  50.95 


inches  of  Hg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  C 
degrees  C 
meters/sec 
meters/sec 


meters/sec 

degrees  C 
degrees  C 


Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =1.72  Kg/m^ 

=  M  =  19.57E-6  Pa.S 
=  C  =  1006.00  J/Kg  °K 
=  k  =  27.80E-3  w/rn^  “k 
=  r  =  0.883 

=  Pr  =  0.708 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m^) 

0.06112 

0.7511 

0.00220 

1171.0 

517 . 18 

30.76 

0.07104 

0.8730 

0.00212 

1320.5 

516 . 71 

30.73 

0.09128 

1.1217 

0.00203 

1614.0 

491.45 

29 . 23 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

NU 

h 

q(Kw/m^) 

0.06112 

0.7512 

0.0034 

1808.4 

822 . 57 

48.9  1 

0.07104 

0.8732 

0.0035 

2163.84 

846 . 76 

50.4  1 

0.09128 

-  -- 

1  .  1220 

0.0033 

2621 . 44 

798 . 38 

47.5  1 

E.12 


RUN  NO  R105 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


=  =90.0 

=  Pj  =  29.3 
=  Pj  =  53.61 
=  Tj  =  23.00 
=  T2  =  80.13 
=  aj  =  344.87 
=  Uj  =  448.90 
=  Mj  =  1.302 
=  Mf  =  1.345 
=  U2  =  153.30 
=  M2  =  0.411 
=  Tj^=  90.15 
=  Tj  =  53.50 


inches  of  Hg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  C 
degrees  C 
meters/sec 
meters/sec 


meters/sec 

degrees  C 
degrees  c 


Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =1.83  Kg/m- 

=  U  =  19.73E-6  Pa.s 
=  C  =  1006.34  J/Kg  °K 
=  k  =  28.04E-3  w/m^  °K 
=  r  =  0.883 
=  Pr  =  0.708 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

NU 

h 

0.06112 

0.8690 

0.00214 

1315.6 

603.6 

0.07104 

1.0101 

0.00207 

1483.9 

585 . 7 

0.09128 

1 . 2979 

0.00197 

1813.4 

557.1 

q(Kw/m^ 


40.5 
39 . 33 
37.41 


RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


X(m) 

Re(  10^) 

St 

NU 

h 

q(  Kw/m^) 

0.06112 

0.8692 

0.0033 

2043 

937 . 3 

62.94 

0.07104 

1.0103 

0.0031 

2217 

875 . 2 

58.75 

0.09128 

1 . 2982 

0 . 0030 

2757 

847.0 

56 . 87 

RUN  NO  R106 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


100 
29 . 3 

54.2 
23 

82 . 3 
344.87 
455.50 
1 . 325 

1 . 369 
162 . 2 
0.43 
93.5 
55.2 


inches  of  Mg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  C 
degrees  C 
meters/sec 
meters/sec 


meters/sec 

degrees  C 
degrees  C 


Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =  1.868  Kg/m^ 

=  11  =  19.80E-6  Pa.s 

=  Cp  =  1006.5  J/Kq  °K 

=  k  =  28. 13E-3  W/rn^  °K 

=  r  =  0.883 

=  Pr  =  0.707 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re( 10S 

St 

NU 

h 

q(Kw/m^)  1 

0.06112 

0.9353 

0.00211 

1394 

641.6 

45 . 23 

0.07104 

1.0871 

0.00205 

1572 

622.6 

43.90 

0.09128 

1 . 3968 

0.00195 

1922 

592 . 2 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

St 

Nu 

h 

- ^ j 

q(Kw/m^) 

0.06112 

0.93717 

0.0032 

2120.3 

975 . 87 

68.8 

0.07104 

1.08933 

0.0033 

2541 . 5 

1006 . 36 

71.0 

0.09128 

1 . 39967 

0.0034 

3315. 1 

1021.61 

72.0  j 

RUN  NO  R107 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


=  P^  =  110 
=  Pj  =  29.3 
=  Pj  =  55 . 1 
=  Tj  =  23 
=  T2  =  88.0 
=  aj  =  344.87 
=  Uj  =  4464 
=  M.  =  1.344 
=  M^  =  1.390 
=  Uj  =  173.3 
=  M2  =  0.455 

Tjjj-  101 

=  Tj  =  58.5 


inches  of  Hg  (Gauge) 
inche'"  ot  ’*  , 
inches  ol  Hg 
degrees  c 
degrees  C 
meters/sec 
mete  <-  s/sec 


meters/sec 

degrees  C 
degrees  C 


Dens i ty 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =  1.918  Kg/m^ 

=  Ji  =  19.93E-6  Pa.s 
=  Cp  =  1006.8  J/Kg  °K 

=  k  =  28.37E-3  w/m"^  °K 

=  r  =  0.883 

=  Pr  =  0.707 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(  10^) 

...  .  .  . 

St 

NU 

h 

q(Kw/m^)  1 

0.06112 

1.01754 

0.00207 

1491 

692 . 3 

54.00  1 

0.07104 

1  .  18276 

0.00201 

1682 

671 . 7 

52.40  I 

0.09128 

1 .51972 

0.00191 

2056 

638.9 

49.83  1 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


X(m) 

Re(10^) 

St 

Nu 

h 

2 

q  ( Kw/m^ ) 

0.06112 

1 

.01795 

0.0035 

2519 

1169 . 3 

91 . 2 

0.07104 

1 

.  18323 

0.0033 

2761 

1102.43 

86.0 

1  0.09128 

1 

.520283 

0.0032 

3440 

1069.03 

83.4  I 

RUN  NO  R108 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Adiabatic  wall  Temperature 
Reference  Temperature 


=  P4  =  120 

=  Pj  =  29.3 
=  P2  =  58.4 
=  Tj  =  23 


“  Cli  “  O  4  •  O 

=  Uj  =  470 
=  Mj  =  1.37 
=  Mj  =  1.41 
=  Uj  =  180.7 
=  M,  =  0.472 
=  T,  =  105.2 
=  T;"=  60.6 


inches  of  Hg| Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  C 
degrees  C 
meters/sec 
meters/sec 


meters/sec 

degrees  C 
degrees  C 


Density 

Dynamic  Viscosity 
Specific  Heat 
Thermal  conductivity 
Recovery  Factor 
Prandtl  Number 


=  p  =  1.951  Kg/m^ 

=  li  =  20.02E-6  Pa.s 

=  C.  =  1007  J/Kg  °K 

=  k  =  28.51E-3  w/rn^ 

=  r  =  0.883 

=  Pr  =  0.707 


RESULTS  BASED  ON  THEORETICAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


x(m) 

Re(10^) 

St 

NU 

h _ 

q(  Kw/m^) 

0.06112 

1 

.0763 

0.00205 

1559.7 

1 

1 

59 .81 

0.07104 

1 

.  2510 

0.00199 

1759.3 

706.04 

58.04 

0.09128 

1 

.6074 

0.00189 

2150.0 

671 .52 

55 . 20 

RESULTS  BASED  ON  EXPERIMENTAL  HEAT  FLUX  AND  MEASURED  MACH  NUMBER 


X(m) 

Re(  10®) 

St 

Nu 

h 

q(Kw/m^) 

0.06112 

1.0765 

0.00311 

2366 . 90 

1104 

91 

94 

0.07104 

1 . 2512 

0.00323 

2857.35 

1147 

0.09128 

1.6077 

0.00285 

3239.40 

1012 

83 

E-16 


APPENDIX  F 


A  summary  of  the  Film  cooling  Runs  is  given  in  the  following 
pages . 
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FILM  COOLING  RUNS  FOR  P^  =  60”  Hq 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Prandtl  Number 


P^  =  60.00 
pj  =  29.07 
Pj  =  44.69 
Tj  =  23.00 
Tj  =  62.48 
aj  =  344.87 
Uj  =  416.80 
Mj  =  1.209 
=  1.286 
U2  =  109.54 
M2  =  0.298 
Pr  =  0.7087 


inches  of  Hg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  c 
degrees  C 
meters/sec 
meters/sec 


meters/sec 


1  Run 
i  No 

. 

Pc 

(psia) 

Pc 

P2 

^2 

Uc 

' 

B 

E 

q  ( Kw/mj ) 
Gage  No  2  | 

1 

§ 

1  FC601 

21 . 76 

1 . 78 

1 . 58 

103.4 

27.4 

0.30 

28.1  1 

FC602 

22.06 

1,79 

1 . 58 

103.4 

55 . 1 

0.60 

27.3  1 

FC603 

22 . 56 

1 . 80 

1.58 

103.4 

82.5 

0.91 

24.6  i 

FC604 

23.27 

1.81 

1.58 

103.4 

109.3 

1.22 

24.2 

■  FC605 

25.41 

1.89 

1 . 58 

103.4 

162.0 

1.85 

23.0  I 

'  FC606 

28 . 63 

1 .92 

1.58 

103.4 

212.5 

2.51 

22.5  j 

;  FC607 

33.20 

2.00 

1.58 

103.4 

260.2 

3.20 

21.2  J 

t  FC608 

39 .43 

2  .  10 

1 . 58 

103.4 

304.6 

3.93 

18.0  1 

P  FC609 

43 . 60 

2 . 27 

1 . 58 

103.4 

312.9 

4 . 35 

17.2  1 

!  FC610 

50.16 

2.61 

1.58 

103.4 

313.2 

5.01 

15.0  ^ 

1  FC611 

59 . 90 

3.11 

1 . 58 

103.4 

313 . 2 

5 . 98 

16.0 

1  FC612 

66 . 39 

3.45 

1 . 58 

103.4 

313 . 2 

6 .63 

17.3  !i 

i!  FC613 

72 . 89 

3 . 79 

1 . 58 

103.4 

313 . 2 

7 . 28 

18.3  i 

'  FC614 

79 . 38 

4.13 

1 . 58 

103.4 

313.2 

7 . 93 

19.0 

1  FC615 

85.87 

4 . 46 

1 . 58 

103.4 

313.2 

8.57 

20.0 

1  FC616 

92 . 37 

4 . 80 

1 . 58 

103.4 

313 . 2 

9 . 22 

21.2 

1  FC617 

1 

li 

IT 

.. 

98 . 86 

5 . 14 

1 . 58 

103.4 

313.2 

9.87 

22,1  I 

J 

i 

ii 
• 

J 
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FILM  COOLING  RUNS  FOR  P^  =  80"  Hq 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Prandtl  Number 


P4  =  80.00 
Pj  =  29.21 
P,  =  50.22 

T.  =  24.00 
T2  =  74.38 
a.  =  345.47 
Uj  =  438.20 
M.  =  1.27  2 
M(  =  1.330 

U,  =  139.22 
M2  =  0.373 
Pr  =  0.708 


inches  of  Hg  (Gauge) 
inches  of  Hg 
inches  of  Hg 
degrees  c 
degrees  C 
meters/sec 
meters/sec 


meters/sec 


i  Run  No 

• 

Pc 

(psia) 

P2 

. 

Pc 

^2 

q  (KW/mj) 
Gage  No  2 

j 

;;  FC801 

24 . 78 

1 . 70 

1.20 

138.2 

40.2 

0.34 

i 

50  1 

i  FC802 

26 . 73 

1.70 

2.04 

138.2 

119.6 

1.04 

46  1 

t  FC803 

28.53 

1 . 70 

2.08 

138.2 

156.7 

1 . 39 

44  ? 

J  FC804 

38 . 16 

1 . 70 

2 . 26 

138.2 

265 . 1 

2.55 

40  1 

FC805 

43 . 23 

1 . 70 

2 . 34 

138.2 

297.9 

2 . 97 

34  1 

1  FC806 

46.21 

1 . 70 

2.39 

138.2 

313.6 

3 . 18 

32  p 

K  FC807 

46 . 86 

1 . 70 

2.41 

138.2 

314 . 9 

3.23 

30  f 

8  FC808 

50.11 

1 . 70 

2.58 

138 . 2 

314 . 9 

3 .45 

27  [ 

1;  FC809 

53 . 36 

1 . 70 

2 . 75 

138.2 

314.9 

3 . 68 

28  ! 

'  FC810 

56 . 60 

1 . 70 

2.91 

138.2 

314.9 

3.90 

29  f 

i  FC811 

66 . 34 

1 . 70 

3.41 

138 . 2 

314 . 9 

4 . 56 

30  “ 

j;  FC812 

12. 84 

1 . 70 

3 . 75 

138.2 

314 . 9 

5.02 

32  1 

:  FC813 

1 

79.33 

1 . 70 

4.08 

138.2 

314 . 9 

5.46 

34  j 

i, 
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FILM  COOLING  RUNS  FOR  =  100"  Hq 


Driver  Pressure 
Driven  Pressure 
Pressure  Behind  the  Shock 
Driver/Driven  Temperature 
Temperature  behind  shock 
Sonic  speed 
Measured  shock  speed 
Measured  shock  Mach  No 
Theoretical  shock  Mach  No 
Flow  velocity  behind  shock 
Flow  Mach  No  behind  shock 
Prandtl  Number 


P4 

m2 

ml 

Ml 

Mt 

^2 

Pr 


100 

inches  of  Hg  (Gauge) 

29.3 

inches  of  Hg 

54.2 

inches  of  Hg 

23 

degrees  C 

82.3 

degrees  C 

344 . 87 

meters/sec 

455.50 

1.325 

1.369 

meters/sec 

162.2 

0.43 

0.707 

meters/sec 

1  Run  No 

1 

Pc 

P2 

Pc 

(psia) 

FC1001 

.  .  . . 

26 . 68 

1 . 78 

2.14 

FC1002 

29 . 13 

1 . 78 

2 . 19 

FC1003 

36 . 48 

1.78 

2.35 

FC1004 

41.72 

1.78 

2.43 

FC1005 

48.60 

1 . 78 

2.54 

FC1006 

53 . 36 

1 . 78 

2.74 

FC1007 

59.85 

1 . 78 

3.07 

i  FC1008 

66 . 34 

1 . 78 

3.40 

f  FC1009 

12. 84 

1 . 78 

3.74 

I  FC1010 

79 . 33 

00 

4.07 

^2 

“c 

B 

■ 

q  f 

( Kw/m, )  E 
Gage  No  1 

^  1 

u 

P 

157.5 

48.11 

0.37 

60 

157.5 

130.12 

1.02 

53 

157.5 

236.68 

1.99 

45 

157 . 5 

270.64 

2 . 35 

40 

157.5 

308.77 

2 .81 

39 

157.5 

314 . 63 

3.08 

42 

157.5 

314.63 

3 .45 

44 

157 . 5 

314.63 

3.83 

46 

i! 

157.5 

314 . 63 

4.21 

48 

P 

[j 

157.5 

314.63 

4 . 58 

51 

1 

t 
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Vita 
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